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ANALYSIS AND TESTING OF TWO-DIMENSIONAL 
SLOT NOZZLE EJECTORS WITH VARIABLE AREA 
MIXING SECTIONS 
By 

Gerald B. Gilbert, PMIip G. Hill 
SUMMARY 

Finite difference computer techniques have been used to calculate the 
detailed performance of air to air two dimensional ejectors with symmetric variable 
area mixing sections and co-axial converging primary nozzles. The successful com- 
pletion of this orogram completes a step in the development of a computer 
program to analyze the ejector of the augmentor wing lift augmentation system for STOL 
aircraft. 


The finite difference computer program analyzes two dimensional mixing 
in converging-diverging jets. The analysis of the primary nozzle assumes correct 
expansion of the flow and is suitable for subsonic and slightly supersonic velocity 
levels. The variation of the mixing section channel walls is assumed to be gradual 

so that the static pressure can be assumed uniform on planes perpendicular to the axis. 

2 

An x-i p coordinate system is used in the solution of the momentum and energy equations 
to remove a singularity condition at the wall. Different assumptions for eddy viscosity 
are made for each distinctly different region of the flow based on information available 
in the literature. 

A test program was run to provide two-dimensional ejector test data for 
verification of the computer analysis. Geometry and primary air operating conditions 
similar to a typical augmentor wing ejector were selected for the tests. A primary 
converging nozzle with a discharge geometry of 0. 125" x 8. 0" was supplied with 
600 SCFM of air at about 35 psia and 180°F. This nozzle was combined with two mixing 
section geometries with throat sizes of 1. 25" x 8.0" and 1. 875" x 8. 0" and was tested 
at a total of 11 operating points. Secondary flow was varied by adding three steps of 
increased restriction to the ejector discharge. For each test mass flow rate, wall 
static pressures and several velocity traverses were recorded for comparison with 
analytical results. 
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The comparisons of wall static pressures, centerline velocity, centerline 
temperature, and velocity profiles between experimental and analytical results at the 
same flow rate were generally very good. The computer program presented in this 
report accurately predicts the performance of the simple two-dimensional ejectors and 
thereby successfully completes the objectives of this program. 
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Section 1 


INTRODUCTION 


1. 1 Background 

The augmentor wing concept under investigation by NASA for STOL aircraft 
lift augmentation is powered by an air to air ejector. The wing boundary layer is drawn 
into the deflected double flap augmentor channel at the trailing edge of the wing and is 
pressurized by a high velocity slot jet which is oriented at an angle to the augmentor 
channel. To predict the performance and to optimize the design of the complete 
augmentor wing, an analytical method is needed to predict the performance of the air 
ejector which powers the augmentor flap section. 

Under contract NAS2-5845 a computer analysis was developed for single 
nozzle axisymmetric ejectors with variable area mixing sections using integral 
techniques^. The ejectors of primary interest in that program and earlier programs 
were high entrainment devices using small amounts of supersonic primary flow to pump 
large amounts of low pressure secondary flow. Good agreement was achieved between 
analytical and experimental results. 

The integral analytical techniques used to analyze the axisymmetric ejector 
configurations are also valid for the analysis of two dimensional ejectors. However, 
the augmentor wing configuration may include asymmetric geometries, inlet flow 
distortions, wall slots, and primary nozzles that are qt large angles to the axis of the 
augmentor mixing section. The integral techniques are not easily adaptable to these more 
complex flows. Finite difference techniques can be used to analyze these more complex 
flow geometries at the expense of increased computer time. 



1. 2 Objectives of Program 


The specific objectives of this investigation are the following: 

(1) to develop a finite difference computer program for the analysis 
of two-dimensional, air ejectors with symmetric variable area 
mixing sections and with co-axial converging primary nozzles. 

(2) to obtain test results with two-dimensional ejector configurations 
so that the analytical methods can be checked. 

By modifying the present analysis additional complicating features of the 
actual augmentor wing ejector may be incorporated into the computer program until the 
complete augmentor wing ejector can be successfully analyzed. 


r 
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Section 2 
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NOMENCLATURE 

(continued) 



P 

ro 
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q 


Nozzle pressure 

Time average static pressure 

Turbulent Prandtl number, — — 
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Turbulent heat transfer, (pv) ' T' 

Gas constant 

Atmospheric temperature 

Time average temperature 

Instantaneous fluctuating temperature 

Jet temperature at the nozzle exit plane 

Flow reference temperature 

Nozzle temperature 

Wall reference temperature 

Time averaged velocity in x- direction 

Instantaneous fluctuating x component of velocity 

Jet centerline velocity at the nozzle exit plane 

til 

Unknown velocity at the n grid point 


Dimensionless velocity in x-direction, 
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Friction velocity, * — 

Time averaged flow velocity in y- direction 
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NOMENCLATURE 

(continued) 
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AX 
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J w 

y+ 
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i b 

Y s 

ip* 

P 

p o 

p* 
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Instantaneous fluctuating y- component of velocity 
Mixing section total flow rate 
Nozzle flow rate 
Secondary flow rate 

Space co-ordinate in the axial direction 
Dimensionless space co-ordinate in the axial direction, 


u x 



Step size in x-direction 

Space co- ordinate perpendicular to axial direction 
Dimensionless space co- ordinate perpendicular to axial direction. 
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Duct half width or duct radius 

y u* 

Dimensionless wall co-ordinate 

Constant, unity for axisymmetric flow and zero for two- 
dimensional flow 

C 

Ratio of specific heat, — 

Cp 

Transformed co-ordinate defined by equation 8 
Regular stream coordinate 

.2 

Dimensionless^ co-ordinate ip *% — — for two-dimensional flow 

^o*o 


Time averaged fluid density 

Fluid density evaluated at a reference temperature T q 
D imensionless fluid density 
Time averaged absolute viscosity 

Absolute viscosity evaluated at a reference temperature T q 

s 

Dimensionless absolute viscosity, - — 

M o 


Mean average shear stress 
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NOMENCLATURE 

(continued) 


Turbulent shear stress, (pv)' u' 
Local wall shear stress 
Eddy viscosity 
Eddy conductivity 


T - T 


Dimensionless temperature 


T - T 
wr o 


Kinematic viscosity at local temperature 

Reference kinematic viscosity evaluated at a reference temperature T 

Local wall boundary layer thickness or jet half width 

u o<5 

Dimensionless boundary layer thickness, — 

o 

Mixing length constant 


Mean value of dissipation 



Section 3 

ANALYSIS OF TWO-DIMENSIONAL JET MIXING 


3. 1 Introduction 

This section is concerned with the essential physical features of a computa- 
tion model for plane two-dimensional jet mixing in converging- diverging jets. A finite- 
difference computer program has been developed for treating the mixing of two parallel 
and compressible air streams, allowing for at least one of them to be supersonic. In 
all cases, the nozzle expansion is assumed "correct", i. e. nozzle exit plane pressure 
is matched to the ambient pressure at that station. Thus, expansion waves and shocks 
at the nozzle exit plane are assumed to be absent. Even though the correct expansion 
assumption may not be realized in a practical case, the downstream flow field will not 
likely be sensitive to small degrees of over - or under- expansion. The flows considered 
include compound flows of supersonic and subsonic streams; however, no provision is 
made for compound choking which may occur with an appropriate transverse distribution 
of Mach number. Such a condition is amenable to analytical treatment under simplified 
circumstances, but has not been encountered in experimental tests carried out so far. 

This development is restricted to symmetric jet mixing in which the high 
speed jet is located on the axis of the channel and no provision is made for blowing or 
suction along the channel walls. The variation in channel geometry along the axis is 
assumed gradual, so that wall curvature is neglected and, on all planes normal to the 
axis, the pressure is assumed uniform. 

In most calculations performed with this method to date, the velocity distri- 
bution at the nozzle exit plane was assumed to be rectangular, i. e. , the wall boundary 
layer has been assumed to have zero thickness at that point; the initial thickness of the 
jet^ secondary stream shear layer has also been assumed to be zero. This requirement 
is not necessary, however, and in general any initial distribution of velocity in the initial 
plane is permissible, under the assumption that pressure distribution across the plane 
is uniform. 

Although previous work ^ has amply demonstrated that integral methods 
are capable of predicting symmetric jet mixing of compressible flow in jets, the finite 
difference method has been chosen for this problem. The finite difference method has 
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advantages relative to the integral method of much greater flexibility in allowable flow 
inlet conditions, and wall boundary conditions, e. g. , the use of wall jets or wall suction. 
Further the finite difference method offers the considerable advantage of mathematical 
precision in determining the overall consequences of any particular physical hypothesis 
regarding the shear stress distribution. With the integral method, the mathematical 
approximation due to the formation of integrals may contribute uncertainty in flow 
prediction in addition to the uncertainty introduced by a lack of precise physical 
knowledge. Thus, in developing a model to handle a certain class of flows, it is 
advantageous to have a method which is relatively precise mathematically, so that the 
effects of physical uncertainties may be assessed relatively clearly. The finite differ- 
ence method is however, quite costly in its requirement for computer time. Further, 
as experience has shown, considerable care is required in adjusting the computation 
grid such that spacings are appropriately small in the region of the wall, and in any 
part of the flow where velocity gradients are quite large. 


3. 2 Basic Conservation Equations 


In 


stream- wise coordinates, the momentum and energy equations 


( 2 ) 


for the 


plane two-dimensional flow are: 


u 9_u 
8 x 


1 dp u 8 r 
? d x d ip s 


( 1 ) 


u8(C p T) 
9 x 


u 

~T 


d x 


u 8 q 

8 4> s P 


$ = 



■) - (P v)' u' 




( 2 ) 

( 3 ) 


in which u is the velocity component in the x or principal flow direction, p is the static 
pressure, p the density and T is the temperature of the fluid. Using the eddy 
viscosity assumption, the mean average shear stress and heat transfer are defined by: 

T = M f-; - <PV)» u' = (P + pe ) fy (4) 


10 



( 5 ) 


i " S T5 - e P «’ v >" r ' * (k + > If- 

in which e is the kinematic eddy viscosity. 

In developing the finite difference solution to this problem, the stream- wise 
coordinate system was attractive, not only in terms of the simplicity of the governing 
equations but also for possible development as a design procedure, in which the flow 
field pressure distribution could be specified and the required wall geometry determined, 
non-interatively, once the solution is obtained in stream coordinates. However, the 
difficulty with the stream wise coordinate is that it introduces a singularity in the 
governing equations in the vicinity of the wall. Given the definition of the stream 
function, 

9* s __ 

Ty = PU (6) 

it can be seen that the gradient 


becomes undefined at the wall where the value of u approaches zero. The singularity 
can be removed as Denny^ ^has shown by using the transformation 

dip 2 - - dip _ P u , 9 u _ p u a u 
9y ^’9y 2 ip * 3y 2i p 9 ip 


instead of conventional stream function definition 

9 u 

in which case the limiting value of the gradient ^ ^ is finite and higher derivatives 
also exist. With this transformation then, the equations of motion may be written. 


- 9 u 
U 9 x 


1 d p 

J> d x 


<M +Pe) 2 “ 
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where ^ Is now the transformed quantity according to Denney^. The transformation 
of these equations is shown in Appendix A. 


3.3 Dimensionless Groups 

Before solution of the finite- difference method, these equations are made 
dimensionless by the following steps. 


The velocity u is normalized by dividing by the jet centerline velocity u 


Also a reference Mach number is defined by: 

u 


M. = 
ir 


yjy^o 


(11) 


in which T q is a reference temperature and y iB the specific heat ratio. A dimension- 
less temperature parameter is defined by: 

T - T. 

9 = 


T - T 
wr o 

in which T wr is a second arbitrary reference temperature. 


( 12 ) 


The fluid properties variables are made dimensionless by defining: 
k 


k* = 


°*p ’ 


n c 

o po 


IS. 
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ro 


po 


(13) 
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in which k Q , Cp Q , P Q > anc * P 0 are fluid properties at reference values of pressure and 

temperature and m • = P v . 

o o o 

In the program the reference values of temperature are 

T - 520°R 
o 

v . T ■= 560°R " ' r ' ” 1 ' " 

wr 

and the reference fluid properties are evaluated at 520°R and 2115 psf. 

The coordinate variables are transformed to: 
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u x 
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(14) 
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V o o 

Then in dimensionless form the equations of motion become: 


(15) 
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L l p* ' { 2 tp* a ip *’ 


(17) 


in which 


C L = 


(y - 1) M. 


ir 


u 


wr 


C (T -T ) 
po' wr o' 


(18) 


The turbulent Prandtl number P ^ is taken to be 0.9. Neglecting the depend- 
ance of the specific heat on temperature, C* = 1. 0. The derivative of the dimensionless 
equations of motion is shown in Appendix A. 
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3.4 Evaluation of the Eddy Viscosity 


In general, the eddy viscosity is evaluated by 


e 


i 2 

m 9 y 


(19) 


in which f is the mixing length. In two-dimensional jet mixing, values of mixing 
length are not well known especially for the region in which the shear zone extends 
from wall to wall. In various zones of the flow, the mixing lengths have been 
evaluated as follows: 


In the shear layer adjacent to the potential core zone of the primary jet 
the mixing length is evaluated from 

f m = 0. 08 6 (20) 

in which 6 is the shear layer width (including the zone between }% and 99% of the 
total velocity difference between primary and secondary streams). 

For the "fully- rounded" portion of the jet flowing coaxially with a 
secondary potential stream, the mixing length has been calculated from 

f m = 0. 108 5 (21) 

in which 6 is the half-width of the jet, evaluated from centerline to the point at 
which the difference between local and secondary velocity is only 1% of the difference 
between centerline and secondary velocity. 

In the wall boundary layer, the mixing length has been evaluated from 
the lesser of: 


f = 0. 096 (outer part) 


or, using the Van Driest approximation, 

(y+/ 26) 


i =0.41 
m 


1 - e 


) 


y (inner part) 


( 22 ) 


(23) 
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in which 



For the region downstream of the point where the jet spreads to inter- 
sect the edge of the boundary layer the mixing length is evaluated, as a first 
approximation only, from 

- y w kl4-0.08 ( ~¥~ ) -0.06 (-f) 4 ] (25) 

L w w j 

which is due to Nikuradse and is cited by Schlichting^ for fully developed flow in 
round tubes. Near the wall (y = y ) the mixing length is evaluated by the Van 
Driest approximation cited earlier, provided the local mixing length so calculated is 
less than that given by the Nikuradse formula. 

3. 5 Boundary Conditions 

With prescribed wall geometry the boundary conditions at the outer wall 

are: 


y = y w ( x ) 

ip* = const. 


d_6 

dip* 


0 


u = 0 


Along the channel axis of symmetry the boundary conditions are: 

y = o 

ip* = 0 
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3.6 Finite Difference Procedure 


By the finite-difference technique, the derivatives in the differential 
equations of motion are replaced by differences either along a streamline between 
two neighboring points X and X + AX or normal to it between two neighboring 
points ip* and <P* + A if>*. 

If one takes the velocity field at plane X as completely known then the 
velocity field at X + AX may be solved, using the implicit method, from the 
finite-difference form of the momentum equation which is of the form 


A n-1 U 2,n + B n-1 U 2, n + 1 + C n-l U 2,n-l D n-1 

th 


(26) 


in which u 9 is the unknown velocity at the n grid point on plane X + AX and 
U) n 

A , B . C , D . are coefficients containing the mean pressure gradient 
n-1, n-1, n-1, n 1 

between X and X + AX and the velocity and shear stress distributions at plane X. 


(5) 

As shown in the derivation in reference' and Appendix B, the co- 
efficients in the finite difference form of the momentum equation are evaluated from: 
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Y8 = 


Y9 = 


u 


hJL 


2f 


u 


liJ 1 


2 tp> 


n 


S . + S 
n+ I n_ 

Afy SI 


S + S. 


n-1 


( A p o SI > 


(31) 

(32) 


16 



SI = A ■» A ip 2 


(33) 


A V *V**n-l’ A *1 ■**«+!-*! 
S “ ( 2^* ^ ^ U 


(34) 

(35) 


,(5) 


In a similar way' , the energy equation can be written in the finite difference form: 

I 

A . 0 + B 0 . + C , 0 O , = D . (36) 

n-1 2,n n-1 2,n+l n-1 2,n-l n-1 
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(44) 
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R 1 ~ A ^ 2 (^P 2 + ) 

and 

A^g 

R2 = a^ <&r 2 + A^) 


(45) 


The relationship between the x-$ coordinates, and the physical plane in 
finite difference form, for any n, becomes, 
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Finally the property relation becomes: 
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(47) 


For a set of N ip-lines and known boundary conditions, Equations (26) 
and (36) each provide a set of N-2 conditions to solve for the unknown velocities and 
temperatures. Each set of equations can be solved simultaneously if the pressure 
gradient is known or assumed. For calculation of flow between fixed channel walls, 
the pressure gradient is assumed and the velocities determined; then the location of 
the outer boundary is calculated from successive use of equation (46) across all N 
grid lines. If the calculated value of the outer boundary location does not agree 
satisfactorily with the actual wall geometry, a new value of the pressure gradient 
is chosen. 

Since each set of equations can be represented by a tridiagonal matrix 
of coefficients, the Thomas Algorithm^ is employed for speedy solution as shown 
in Appendix C which describes the solution procedure. 

The structure of the computer program is given in Appendix D. 
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Section 4 


TEST PROGRAM 

A two-dimensional experimental rig was designed, fabricated, and installed 
in our laboratory. The purpose of the experimental work was to obtain test data for 
verification and adjustment of the computer analysis. The experimental program is 
described in this section. 

4. 1 Experimental Apparatus 

4. 1. 1 Two-Dimensional Ejector 

The two-dimensional ejector consisted of a slot type primary nozzle and a 
two-dimensional mixing section. The arrangement of the ejector system is shown on 
Figure 1. 

A picture of the primary nozzle is shown on Figure 2. The discharge slot 
is 0. 1215" + . 0005" by 8. 00" with rounded corners. The side walls are quarter inch 
carbon steel and four internal supports are included to prevent widening of the discharge 
slot when the nozzle is pressurized. Dial indicator measurements show that the slot 
opened up by about 0. 0008 inches in the center of the nozzle, about . 0004" at the 
quarter width location and zero near the ends of the slot. This is equivalent to an 
increase in nozzle slot area of 0. 33% when pressurized. Stagnation pressure measure- 
ments were made with a kiel probe from side to side in the nozzle discharge and were 
found to be uniform across the 8" width of the slot. The primary nozzle is positioned 
in the mixing section (see Figure 1 and Figure 3) so that the primary flow is discharged 
along the centerline of the straight symmetrical mixing section. 

The mixing section as shown on Figure 1 consists of a rectangular variable 
area channel formed by two identically contoured aluminum plates and two flat side 
plates. The pictures in Figures 3 and 4 show two views of the mixing section. The 
two contoured plates can be positioned in two symmetrical locations about the center- 
line to form the two channels tested (throat heights of 1. 25" and 1. 875"). The width 
of the mixing section is 8. 00" for the full length. The variation of channel height with 
distance from the nozzle discharge is given on Table 1 for the 1. 875 throat mixing 
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section. The geometry for the 1. 25" throat height is obtained by subtracting 0. 312" 
from each y value. Three plexiglass windows are installed along each side of the 
mixing section so the tufts of wool mounted inside can be observed for indications of 
flow separations and unsteadiness. 

The screened mixing section inlet is shown on Figure 5. Initial tests with- 
out the extended inlet showed that highly swirling corner vortices were formed in the 
four corners of the bellmouth and extended into the test section. The extended inlet 
eliminated the corner vortices and improved the stability of the ejector flow and static 
pressures. The extended inlet shown on Figure 5 was used for all ejector tests. 

4.1.2 Facilities for Ejector Tests 

The schematic of the ejector test facilities on Figure 6 shows the three 
required subsystems needed for operation, control and measurement of the ejector: 


* Primary Flow System 
« Mixed Flow System 
« Boundary Layer Suction System 

The primary air flow is supplied by a 900 SCFM oil free screw compressor 
at 100 psig and an equilibrium operating temperature between 180°F and 240°F. The 
primary air flow rate and pressure are controlled by a manual pressure regulator and 
bleed valve. The mass flow is measured by a standard 3 inch Danial orifice system. 
The air flow is delivered to the primary nozzle through a flexible hose. 


The mixed flow system consists of a plenum chamber, an 8" orifice 
system and a throttle valve. Four different operating flow rates are achieved by the 
following equipment combinations. 


1. Maximum Flow Rate - 

2. First Reduced Flow Rate - 

3. Second Reduced Flow Rate - 

4. Lowest Flow Rate - 


Mixed flow discharges directly into 
laboratory from mixing section. 

The plenum is connected to the mixing 
section discharge. 

The orifice is connected to the plenum. 
The throttle valve is partially closed. 
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Orifice flow rates are obtained only for the two lowest flow rate conditions. Figure 7 
and 8 show most of the experimental ejector installation. The large rectangular box 
connected to the mixing section by the large black flexible hose is the main plenum. 

The 8" orifice is not visible in the picture. 

The suction system removes the boundary layer flow from each of the four 
corners of the mixing section to prevent wall boundary layer separation in the ejector. 
The pictures in Figures 7 and 8 show three 3/4 inch tubes connected to each corner of 
the mixing section. These 12 tubes collect the boundary layer flow from the corner 
suction slots which are 0. 060 inches wide and are machined into the sides of the 
contoured plates (See figures 9 and 10). The four tubes at one X location are connected 
to a single large tube under the mounting table. The three large tubes are each 
connected to a large tank plenum through a separate throttle valve. A Roots blower 
draws the air through the suction system and through a three inch orifice system. The 
suction system is capable of removing about 1% to 2% of the mixing section flow rate. 
During the operation of the ejector rig, the boundary layer suction system was 
necessary to prevent flow separation in the mixing section diffuser. The presence of 
separation was easily observed from the violently flopping tufts, the large fluctuation 
in wall static pressures and audible pulsations. The operation of the suction system 
drastically reduced these symptoms. 

The ejector system was operated by starting the primary air flow at low 
pressure and flow rate. The suction was turned on and then the primary pressure 
was increased to the desired test conditions. The large mixing section (1. 875" 
throat height) was operated at 21 psig without separation in the mixing section. The 
small mixing section (1. 25" throat height) could not be operated over 20 psig without 
separation for the high flow condition. The tests with the small mixing were therefore 
run at 17 psig. 

4. 2 Instrumentation and Data Reduction 
4. 2. 1 Instrumentation 


The following instrumentation was included on the test rig. 
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Primary Flow System 


Flow Rate - Standard 3" orifice system 

Nozzle Pressure - Pressure gage accurate to + .25 psig 

Nozzle Temperature - Thermocouple with digital readout 


Mixed Flow System 


Flow Rate - 8" orifice system for two lowest flow rate conditions 

Static Pressures - Wall static pressures down the center of the mixing 
section and some at other locations (see Figures 9 
and 10). Manometers were used for measurement. 
Traverse Data - Stagnation pressure and temperature profiles were 

measured at up to 9 axial locations using a kiel temper- 
ature probe, a pressure transducer and direct digital 
readout, and a temperature direct digital readout 
(see Figure 8). 


Suction Flow System 

Flow Rate - 3" orifice system 

Suction Pressure- a mercury manometer 


4.2.2 Data Reduction Procedures 

Three types of data reduction calculations were needed in this program: 


• Standard orifice calculations 

• velocity profile calculations 

• integration of velocity profiles to calculate flow rate 

The orifice calculations were carried out using standard orifice equations and ASME 
orifice coefficients. The velocity profiles were calculated from the well known 
compressible flow relationships between Mach number and the ratio of stagnation 
pressure to static pressure that can be found in most fluid mechanics text books. The 
local velocity is calculated from the Mach number and the local speed of sound which 
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is dependent on the local static temperature. The static temperature is calculated from 
the measured stagnation temperature profiles and the compressible flow relation be- 
tween temperature ratio and Mach number. 

To calculate an integrated mass flow rate for each traverse location a time 
sharing data reduction computer program was written to integrate the product of local 
velocity and local density over a two-dimensional section of unit width. The program 
also calculated the "mass-momentum" stagnation pressure at each traverse section 
using the equations presented on page 52 and 53 of reference 6. The mass-momentum 
method determines the flow conditions for a uniform velocity profile which has the same 
integrated values of mass flow rate, momentum, and energy as the non-uniform velocity 
profile actually present. 

4.2.3 Experimental Uncertainty 

Orifice Calculations 


The techniques presented in reference 7 were applied to the primary flow 


orifice calculations and the mixed flow orifice calculations. 

The following uncertainty 

results were obtained: 



Orifice 

Nozzle Pressure 

Uncertainty 


psig 


Primary 

17. 0 and 21 

+ 0. 8% 

Mixed 

slightly above 



atmospheric 

+ 1. 3% 

Static Pressures 




Uncertainty in the wall static pressures mainly occurs because of un- 
steadiness in the manometer liquid columns caused by unsteadiness in the flow. The 
lowest flow rate condition which had the most system resistance downstream of the 
mixing section had a wall static pressure unsteadiness of about +3/8 inches of water. 
The amount of unsteadiness increased as the flow rate was increased by removing 
system resistance. For the unrestricted maximum flow rate condition the wall static 
pressure unsteadiness was + 2.0 inches of water. These values are also a measure 
of the uncertainty. 
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Integrated Mass Flow Hate 

The mass flow rate calculated by integrating the results of the stagnation 
pressure and temperature traverses is influenced by many items and is therefore 
very difficult to estimate. The following items all contribute to the uncertainty in 
integrated mass flow rate: 

1. unsteady wall static pressures 

2. unsteady traverse stagnation pressures 

3. instrument accuracy of the pressure transducer and digital readout 

4. inaccuracies due to the effect of steep velocity gradients on sensed 
pressure 

5. inaccuracies due to probe effect near the mixing section walls 

6. inaccuracy in probe position 

7. assumptions and inaccuracies associated with the data reduction 
computer program 

8. data recording errors or computer data input errors 

9. errors caused by loose connections in the pneumatic sensing tube 
between the probe and the transducer 

10. Non- two- dimensional flow distribution across the width of the 8 inch 
mixing section. 

All of these effects could combine to give both a + uncertainty band and a fixed error 
shift. 


One measure of the uncertainty due to these effects is obtained from the 
limits of individual integrated mass flows for each test run. These values are listed 
on Table 2 for all of the test runs with traverse data. The results presented on 
Table 2 show an average variation of + 3. 6% and -2. 8% or a total spread of 6. 4%. 
These values only include the effect of variable uncertainty and exclude the uncertainty 
due to probe errors in steep gradients and near walls and integration assumptions. 
Both of the excluded errors probably cause the intergrated mass flows to be too large 
because the probe tends to measure too high near the wall and the integration program 
neglects wall boundary layers. 
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From the above discussion it is concluded that the average integrated mass 
flow rates may have a fixed error of +1% to 2% and an uncertainty of about +3% to +4%. 

4.3 Test Results 

A total of eleven ejector tests were carried out on two mixing section con- 
figurations (1. 25" and 1. 875" throat height). The data presented in this report falls 
i nto the following categories: 

Test Conditions and Mass Flows 
Static Pressures 

Centerline Velocities and Temperatures 
Velocity Profiles 
Temperature Profiles 
Eddy viscosity Sensitivity 
Flow Rate Sensitivity 

Table 3 shows which figures and tables show the data for each test run. Most of the 
figures and tables present both test data and comparative analytical results. The 
comparisons will be discussed in section 5. 0. 
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Section 5 


COMPARISON OF ANALYTICAL AND EXPERIMENTAL RESULTS 
5. 1 Test Conditions and Mass Flows 


Table 4 presents a tabulation of the measured nozzle conditions, the inte- 
grated mass flow rate from the measured pressure and temperature profiles, and the 
integrated "mass momentum" stagnation pressure. 

The nozzle mass flow rate was calculated from standard orifice readings 
which were shown in section 4. 2. 3 to have an uncertainty of about +0. 8%. Using the 
orifice flow rate, the nozzle pressure, the nozzle temperature, and the nozzle discharge 
area, a nozzle discharge coefficient (C^) was calculated for each test run. These values 
all fall within a range of +0. 007 and -0. 0085 around an average of 0. 973 which is con- 
sistent with the calculated uncertainty. If there were no error in the nozzle calculations 

all of the CL, values would be identical. From these results it is safe to assume that 
N 

the listed nozzle flow rates are accurate to at least +1%. 

The tabulated mixing section flow rates were calculated as described in 
section 4. 2. 2 by integrating the measured pressure and temperature profiles. As 
described in section 4. 2. 3, these results probably have a fixed error of between +1% 
and + 2 % and an uncertainty of between +3% and +4%. Table 5 presents a comparison 
between three separate mass flow determinations: 

• integrated from traverse data 

• measured by orifice 

• computer mass flow giving the best wall static pressure 
comparison 

Only 4 of the tests could be measured with the large orifice, but all of these four tests 
agree with the computer mass flow within _+0. 9% as shown on table 5. Section 4. 2. 3 
shows that the expected uncertainty in orifice mass flow is about +1. 3% making it much 
more accurate than the integrated traverse values. The wall static pressures are in 
fact a function of the average mass flow represented by the orifice value rather than a 
local velocity profile down the center of the two-dimensional mixing section. This is 
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true because the mixing section flow patterns can not support a side- to- side pressure 
gradient along the 8 inch width of the mixing section which was verified by test measure- 
ments. Therefore it is concluded that the measured orifice mass flows and the computer 
mass flow for best match of wall static pressures are the correct mass flow values. 

The integrated mass flows are in error and in some cases inconsistent. Table 5 shows 
that the integrated mass flow values spread over a range of -2. 9% to +6. 4% around the 
computer determined value. Figures 11 and 12 show all of the mass flow values on 
Table 5 plotted versus the mixing section throat static pressure. Figure 11 for Runs 
1-5 shows the good agreement between computer analytical mass flows and orifice mass 
flows and shows the wide scatter of integrated traverse mass flows. Figure 12 for Runs 
6- 10 again shows good agreement between analytical and orifice values and this time 
shows a consistent trend of integrated traverse mass flows which are now offset by 
about +3. 2% on a line parallel to the other more accurate mass flow values. 

The "mass- momentum" stagnation pressure listed on table 4 suffers from 
the same inaccuracies as the integrated mass flow rate discussed above. The plotting 
of mass-momentum stagnation pressure versus mass flow will therefore show some 
discrepancies. 

5. 2 Mixing Section Wall Static Pressure Variation 

The wall static pressure distributions are shown on Figures 13 and 14 and 
Table 6 as specified on Table 3. Runs 4, 8, and 11 on Table 6 were extra tests for 
which no analytical solutions were obtained. Test Run 11 was a repeat of test Run 9 
and gives results that are essentially the same. 

Figures 13 and 14 show there is a good comparison between experimental 
wall static pressures (shown as data points) and the analytical static pressures (solid 
lines) at essentially the same mass flow (see discussion in section 5. 1). The analytical 
results have assumed that the mixing length constant in equation 20 is 0. 08 and in 
equation 21 is 0. 108. These values influence the mixing process through the eddy 
viscosity. The influence on wall pressures is relatively minor as will be discussed in 
section 5. 5 where these values are varied over a reasonable range. The comparison 
between test and analytical values is generally excellent. Both the data and analytical 
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results show changes in shape at points where the geometry changes. The two areas 
where some disagreement occurs is in the entrance region and in the last half of the 
diffuser. 


The difference in the bellmouth section occurs because the analytical 
program calculates a centerline static pressure and assumes the static pressure con- 
stant at each x distance from the nozzle discharge whereas the experimental data are 
wall static pressures and can be influenced by curving streamlines. At x = 0 the bell- 
mouth walls still have a significant curvature which causes flow streamline curvature 
in this region. The result is a reduced wall static pressure and an elevated centerline 
static pressure. Between 1 and 2 inches downstream of the nozzle discharge the wall 
curvature is reduced to very small values and the data and analytical results agree very 
closely. 


The second area where minor differences occur is in the last half of the 
diffuser for the higher flow rate test runs. The reason for this difference could be an 
underestimation of the pressure losses due to wall friction, mixing, and diffusion. 
Substantiation of this can be seen by comparing the slope of the pressure data to the 
analytical results in the constant area throat section between 8 and 11 inches. For the 
low flow rate Runs 2, 6 and 7 where the slopes are essentially equal, the test and 
analytical diffuser wall pressures are almost identical. For the other runs the test 
data slope between 8 and 11 inches is always more negative than the analytical results. 
For frictionless uniform flow in a short constant area duct, the static pressures would 
be equal all along the duct. For frictionless non-uniform flow in a short constant area 
duct the static pressure can increase as mixing takes place. For non-uniform flow in a 
constant area duct with friction, the static pressure will tend to decrease along the duct 
and the slope will become less positive or more negative as flow rate (and therefore 
losses) increases. From these observations, it would appear that the flow dependent 
losses for the analytical solution may be underestimated in the constant area and 
diffusing sections. This may be the cause of the difference between the test and 
analytical wall static pressures in the diffuser section. 
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5. 3 Centerline Velocity and Temperature Variations 

Figures 15 and 16 present the variation of maximum velocity and maximum 
temperature as a function of distance from the nozzle discharge. The temperature 
comparison is generally good for all test runs. The velocity comparison is also good. 
However the experimental maximum velocities tend to be higher than the analytical 
values in the first 4 inches downstream of the nozzle discharge. In the throat section 
and diffuser, the experimental values tend to be lower than the analytical values. In 
general the comparisons are very good. Differences may occur due to the eddy viscosity 
and mixing length distributions assumed (see section 3.4) or due to measurement in- 
accuracies. , 

5. 4 Velocity Profiles and Temperature Profiles 

A total of 45 sets of traverse measurements were taken during the experi- 
mental test program. Table 3 shows the figure numbers that present the comparison 
of the test data and analytical results for each test run. These results are presented 
on Figures 17 through 26. 

In general the comparison of profile shape and velocity magnitude is very 
good between the analytical and experimental profiles. The comparisons for Runs 6 
through 10 (Figures 21-24) match very closely. The only differences that are notice- 
able are that the experimental velocity profiles within 5. 0 inches of the nozzle discharge 
are off center by about O'. 025" and slightly higher in maximum velocity than the corres- 
ponding analytical velocities. The nonsymmetry has disappeared for all traverses at 
distances greater than 5 inches. The good match of velocity profiles for Runs 6 through 
10 goes along with the good comparison of static pressures and the consistent trend in 
integrated traverse mass flow rate discussed previously. 

The comparison of experimental and analytical velocity and temperatures is 
not as good for Runs 1 through 5 as it was for Runs 6 through 10. The comparisons are 
also not as consistent from run to run which also coincides with some of the static 
pressure and mass flow differences noted previously for these runs. The following 
observations apply only to Runs 1 through 5. 
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1. The experimental jet is off center by about 0. 057" but the 
non- symmetry has disappeared for profiles at distances 
of greater than 5.0". 

2. For x of 3. 0" or less the peak experimental velocities 
are greater than the analytical values for Run 3 and 
Run 2 and are slightly less for Runs 1 and 5. 

3. The spread width of the velocity profiles compares very 
well at distances from the nozzle of 7. 0 inches or less. 

For distances between 7 inches and 16 inches, the 
experimental profiles tend to spread faster and have a 
flatter profile. 

4. The experimental temperature profiles in Figure 25 are 
spread significantly more than the analytical values at 

x = 3. 0" and x - 10. 5", the only two profiles plotted. 

5. The comparisons for Run 1 are better than for the other 
runs for the 1. 25" throat mixing section. 

Both sets of data (for the 1. 25" and 1. 875" throat height) were calculated 
using the same eddy viscosity assumptions for mixing (0.08 for eq. 20, 0. 108 for 
eq. 21). The test Runs 6 through 10 have lower average throat Mach numbers (. 39 to 
.52), slightly higher primary nozzle velocities, higher wall static pressures, and 
larger mixing section dimensions. The eddy viscosity assumptions may be more suit- 
able for these operating conditions than for those of test Runs 1 through 5. In any 
event, the agreement between experimental and analytical results is better for the 
Runs 6 through 10. 

5. 5 Sensitivity of Computer Analysis 

The sensitivity of the computer analysis to changes in eddy viscosity and 
flow rate were investigated to obtain a measure of the amount of performance change 
that can result from small changes in assumed values. 
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5.5.1 Eddy Viscosity 

The results for the eddy viscosity changes are shown on Figures 27, 28, and 
29. The eddy viscosity is directly proportional to the square of the mixing length 
according to equation 19. The changes in mixing length were confined to the mixing 
region prior to the point where the jet mixing reaches the developing wall boundary 
layer. In this region the mixing length is defined by equations 20 and 21 as a constant 
times a mixing zone dimension (see section 3.4) 

Equation 20 is used to calculate the mixing length in the region close to the 
nozzle discharge where the primary jet still has a flat potential core (probably confined 
to the first 0. 5" to 1. 0" of mixing). Most of the calculations have been carried out 
using a constant of 0. 08 in equation 20. For the results presented in this section the 
comparative runs were made with the constant equal to 0. 094 which gives about a 38% 
increase in eddy viscosity in this small region. 

Equation 21 is used to calculate the mixing length in the region where the 
primary jet is "fully rounded" but has not intersected with the wall boundary layer. 

This region extends for about 4" to 6" into the mixing section for the 1. 25" throat con- 
figuration and extends for about 6" to 8" for the 1. 875" throat configuration. Most of 
the calculations have been carried out using a constant of 0. 108 in equation 21. For the 
results presented in this section, the comparative runs were made with a constant equal 
to 0. 120 which gives about a 23% increase in eddy viscosity. 

The velocity and temperature results shown on Figures 28 and 29 for Runs 
3 and 6 show that the amount of mixing increases with eddy viscosity. This results in 
reduced centerline velocities, increased velocities near the walls and increased wall 
static pressures (see Figure 27). All of the changes are small. 

The effect of mixing length changes in the rest of the mixing section as 
defined by equation 25 was not investigated but it is expected that the results would be 
similar. Section 3.4 points out that equation 25 was obtained by Nikuradse for fully 
developed flow in round tubes and should be considered to give only approximate results. 
Changes in this equation could provide a better match of static pressures for some of 
the high flow test runs as discussed in Section 5.2. 
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5. 5. 2 Flow Rate 


Figure 27 shows the effect on wall static pressures of a 2.2% change in 
total mass flow for Runs 3 and 6. The wall pressure decrease as flow rate is increased 
is about double for Run 3 as compared to Run 6. The reason for this is that the average 
Mach number for Rim 3 (1.25" throat) is larger than for Run 6 (1.875" throat) even 
though the Run 6 mass flow is larger. Figure 30 shows the influence of throat Mach 
number on throat static pressure level. The local slope of this line indicates the rate 
of change of throat pressure with Mach number. Rim 6 happens to be the lowest Mach 
number test run and Run 3 has one of the largest Mach numbers. A 'comparison of the 
local slopes for Run 3 and Run 6 on Figure 30 gives results consistent with Figure 27. 
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Section 6 


CONCLUSIONS 

(1) The finite difference computer analysis developed to analyze 
two-dimensional co-axial slot ejectors with variable area 
mixing sections predicts the performance of the experimental 
configurations tested under this program very closely. 

(2) The analytical and experimental results compared are at 
essentially the same flow rate within the accuracy of our 
measurements. The correct mixing section mass flow 
rates for each test are best represented by the orifice 
measured values and the computer analytical mass flow 
for best comparison of measured wall static pressures. 

These two values agree within - 0. 9%. The integrated 
traverse mass flows are less accurate and range between 
-2.9% and 6. 4% of the other values. 

(3) The experimental and analytical wall static pressure distributions 
agree within 1 or 2 inches of water over most of the mixing 
section for most of the test runs. 

(4) The experimental and analytical velocity profiles compare 
very well in both velocity level and amount of jet spread due 
to mixing. 
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Appendix A 


BASIC EQUATIONS OF MOTION 

The momentum and energy equations as shown in equations 1 and 2 in the 

2 (3) 

main text can be transformed to the xh/> coordinates according to Denny v ’ by the 
following steps. 

Momentum Equation : 

The stream function transformation is defined by: 
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The third term of the momentum equation becomes: 
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The substitution of equation A- 4 into equation 1 of the main text results in 
equation 9 of the main text. 

Energy Equation 

The third term of the energy equation (equation 2) is transformed as follows: 
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The substitution of equation A- 6 into A- 5 completes the transformation of the third term 
of the energy equation as shown in equation A- 7. 
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The fourth term of the energy equation (equation 2 and 3) is transformed by 
substituting equation A- 2 into equation 3 as follows: 
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The substitution of equations (A- 7) and (A- 8) into equation 2 of the main text results in 
equation 10 of the main text. 


Dimensionless Momentum Equation 

The equations 11 through 15 of the main text define the dimensionless groups 
used to non-dimensionalize both the momentum and energy equations. 


The first term of the momentum equation (equation 9) is nonr-dimensional- 
ized as follows: 
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The second term of the momentum equation is non-dimensionalized as 

follows; 
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The third 

term of the momentum equation is 

non-dimensionalized as 



follows: 
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The non-dimensionalized form of the momentum equation (equation 16) is 
obtained by substituting equations A- 9 , A- 1 1 , and A- 12 into equation 9 of the main text 

3 

and eliminating the factor ( u Q /^ 0 ) from each term. 

Dimensionless Energy Equation 

The first term of the energy equation (equation 10) is non-dimensionalized 


as follows: 
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The second term of the energy equation is non-dimensionalized as follows: 
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The third term of the energy equation is non-dimensionalized as follows: 
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The fourth term of the energy equation is non-dimensionalized as follows: 
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Each of the four terms of the energy equation is then divided by the 


(A- 16) 
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(A- 17) 


which results in the following combination of quantities in the second and fourth terms 
of the energy equation: 
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which equals 


o 

C (T -T ) 
po' wr o 

The substitution of equations A- 13, A- 14, A- 15, and A- 16 into equation 10, 
the division by the quantity in (A- 17) and the substitution of into the second and fourth 
terms results in equation 17 of the main text. 
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Appendix B 


FINITE DIFFERENCE EQUATIONS 


This Appendix provides the detailed derivations of the finite difference 
equivalents of the momentum and energy conservation equations, (16) and (17) 
respectively. For convenience the following definitions are introduced: 



These definitions and the assumption that C* = 1. 0 permit the momentum and energy 
equations to be expressed as 


a u _ 
a x 

1 dP 
2p* dX 
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[aul 
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fB- 2\ 

ax 

2p* dX 
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2ip* 

a i ip * 



Before approximating these equations with finite difference relations a 

system of grid lines parallel to the X and ip* axes must be introduced. As illustrated 

in figure B-l , a nodal point coincides with each intersection of these lines. Lines 

parallel to the ip* axis are termed m-lines and those parallel to X axis n- lines. Each 

node is given a double subscript, the first being the number of the m-line passing 

through it, and the second the n-line number. 

m= 1 m=2 n\= 3 



Figure B-l Definition of Grid Lines for Finite Difference Solution 
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The values of the variables on the m=l line are the known initial conditions. The 
conservation equations express for each node on the m=2 line its inter-relation with 
other nodes on the m=2 line and nodes on the m=l line. If m=2 line nodes are only 
related to nodes which lie on the m=l line, the finite difference scheme is termed 
explicit. If an m=2 node is also related to a number of other m=2 nodes, the scheme is 
termed implicit (See figure B-2). 



Figure B-2 

Diagrams of Explicit and Implicit Solutions 

The implicit form of finite difference schemes leads to a series of N simultaneous 
algebraic equations relating the known initial conditions on the m=l line and the unknown 
variables on each of the N nodes on the m=2 line. After solution of these simultaneous 
equations, the variables on the m=3 line are expressed in terms of the known values on 
the m=2 line. Proceeding in this manner, a solution to the complete flow field is 
marched out. Although simpler to program, the explicit scheme shows unstable 
characteristics if the m- lines are widely spaced relative to the n-line spacing. Implicit 
schemes show much more stable characteristics and therefore allow much larger m-line 
spacings, thus reducing computation times. The computer procedure presented in this 
report employs a system of implicit finite difference approximations which are defined 
using the notation described in figure B-3. 


40 





Figure B-3 

Implicit Finite Difference Term Definition 


The velocity at nodes n+1 and n-1 can be expressed in terms of a Taylor series expanded 
about node n, on the same m-line, 
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Neglecting terms of the order (Ai p) and higher, yields 
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The second derivative term in the momentum equation is 
approximated using the following Taylor series expansions. 
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Neglecting terms of the order of — -jp- and higher yields. 
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Similarly, 
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The velocity at a node located at the intersection 
of the downstream m-line and any n-line u 2,n can eX P resse< ^ 
in terms of the following Taylor series , 

(AX) 2 + higher order terms 
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(B-U) 
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Use of the boundary layer equations implies that gradients 

in the X -direction are much smaller than those in the ^-direction. 

Therefore it is permissible to use a simplier approximation of 

the X-direction derivatives. 

Neglecting terms of (AX) 2 and higher yields. 
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This approximation is termed "backward-difference . 
Similarly , 


A — 0 

3 G 1 _ 2 ,n l,n 

3x| n - AX 


(B-12) 


(B-13) 


44 



The only terms in the energy and momentum equations which cannot be approximated 
using the preceeding equations are those containing the pressure gradient . 
Assuming this gradient varies linearly throughout the AX interval yields, ^ 


dP 

dX 



(B- 14) 


Momentum Equation 


Combining equations (B-l), (B-9), (B-12) and (B-14) yields 



This equation can be expressed in the form 


A n-1 U 2,n + B n -1 u 2 ,n+l + C n-1 U 2,n-1 ~ D n-1 (B-16) 

in which the coefficients are defined by equations 27 through 34 of the main text. 

Energy Equation 


Combining equations (B-2), (B-5), (B-10), (B-13) and (B-14) yields 
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This equation can be expressed in the form 
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2, n 
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(B-18) 


in which the coefficients are defined by equations 37 through 45 of the main text. 
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Appendix C 


Solution Procedure 

The calculation procedure starts at the upstream flow boundary, where the 
values of all flow variables must be known or assumed. Specification of the velocity and 
temperature distribution, dimensionless eddy viscosity, duct and nozzle inlet 
dimensions, and working fluid, defines all initial conditions. 

The known initial conditions, m=l line, are related to the unknown 
conditions, m=2 line, by the previously derived equations, and assumed boundary 
conditions. These inter-relations form a set of N-2 simultaneous algebraic equations, 
where N is the number of n- lines, and the equations are shown in Appendix B. The 
resultant matrix of coefficients is tridiagonal in form except for the initial and final rows 
which only contain two terms. Rapid, exact solutions to this type of matrix are ob- 
tained using the Thomas Algorithm, a successive elimination technique, which is 
described in this Appendix. 

The solution for the variables on the m=2 line is iterative, because of the 
presence of the unknown pressure in the momentum equation. The procedure adopted 
was to estimate the pressure gradient, and solve the equations, using the algorithm. 

The equations automatically satisfy conservation of mass, momentum, and energy, but 
only one pressure gradient yields the correct wall geometry. The duct dimension 
corresponding to the estimated pressure gradient was calculated from the m=2 line 
variables. The pressure gradient was then incremented by a small percentage of its 
initial estimated value, and the calculation process repeated for a new duct dimension. 

A third estimate of the pressure gradient was obtained by interpolation between the two 
calculated, and the actual duct dimension. In almost all the calculations performed to 
date, this value has been acceptably close, within 0.001%, to the actual duct dimension. 

If this criterion is not met, a further iteration is applied, and a fourth solution obtained. 

The now known variables on the m=2 line become the new m=l line variables 
and the procedure is repeated for another set of m=2 line variables. Thus a solution to 
the complete flow field is marched out. 



The difference form of the momentum and energy 


equation is: 


Vl X n + B n-1 Vl + C n-1 X n-1 ' D n-1 

^ is either u or 0# If the number of n^lines is N, 
there are N-2 equations of the form (1) and two equations 
expressing the boundary conditions * The first and the last 
equations represent the boundary conditions, which in difference 
form along the axis of symmetry are: 
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(C-3) 


Equations (C-2) and (C-3) can be written in terms of X as follows: 


X 1 = X 2 


(C-4) 


At the duct wall the boundary conditions are 




(C-5) 
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(C-6) 


0 2,N 0 2,N-1 ’ 

Equations (C-5) and (C-6) can be written in terms of X as follows: 

*N = K *N-1 (C-7) 

where K is 0 for the momentum equation and unity for the energy equation. Thus, the 
matrix form of the equation (C— 1) is shown on the following page (Table C-l). 

The second equation is 

C l X l + A l*2 + B 1 X 3 ’ D 1 

| 

Substituting equation (C-4) into this equation yields: 

A X X 2 + b i X 3 = 

where A ' ^ = c 1 + 

The N th -1 equation is 

C n-2 ^N-2 + *N-2 Vl + B N -2 *N = D N -2 

Substituting equation (C-7) into this equation yields: 

I , r •' 

C N-2 V-2 + A N-2 Vl " d n-2 


where a n _ 2 = ^2 + K - B N -2 

Thus the N equations (C-8) can be reduced to the N-2 equations shown on Table C-2. 


(C- 11) 
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(C-9) 


(C-10) 
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ated As Equation C-8 
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The Thomas Algorithm 


Starting with the first equation, X 2 can be 
expressed in terms of X The seconci equation gives X^ in 
terms of X^. Continuing through all the equations until the 
N th -3 equation gives X N _ 2 in terms of X^-^. Combining this 
with the last equation gives X N _^. Working backwards through 
the equations then allows the remaining unknowns to be found. 
This procedure is most easily applied by defining the following: 


I 



n = 2,3 (N-2) (C-14) 

n = 2,3 (N- 2) 

n = 2,3 (N-2) 


Equations (C- 13) then reduce to: 

*0-1 = %-2 and X n = 9 n -l ‘ Q n-1 X n-H n " (N ~ 2) ' (N-3) , 2 

(C-15) 

If the values of W, Q and g are calculated in order of in- 
creasing n using equations (C-14), then equations (C-15) can be 
used to calculate the values of X in order of decreasing X 
starting with X N _^ . To clarify this procedure, the method 
is now used to solve the following four simultaneous equations: 
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hence X 4 = X 5 



(C-17) 



hence X 5 = g 4 (C-18) 

Substituting in equation (C-16) yields X g . Equations (C-17) and (C-18) are special 
forms of equations (C- 15) for N=b and n=4. 
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Appendix D 

COMPUTER PROGRAM 

The computational procedure consists of a main program, which is divided 

into ten sections, and six subroutines. The program Flow Chart is shown on Figure D-l. 

The functions of each section of the main program are as follows: 

Section (1) : Input and Initialization 

(a) Constants which have single, initial value for most applications are defined 
with data statements. 

(b) The parameters which specify the test conditions are inputed from data 
cards. 

(c) Dimensional parts of the data are non-dimensionalized. 

Section (2) : Initial Profiles Generated 

(a) The initial u, 9 , p* , p* , E, Y, and distributions are calculated. 

(b) The shear layer and wall boundary layer thickness are calculated using a 
search technique applied to the m = 1 line velocity profile. 

S ection (3) : Turbulence Model 

(a) The dimensionless eddy viscosity, which will subsequently be used in 
calculating the variables on the m = 2 line, is calculated from m = 1 
velocity profile and one of the turbulence models detailed in the main 
text. 

Section (4): Choice of X-Step 

(a) The distance between the m = 1 and m = 2 lines is chosen. Initially, 
this distance is related to the shear layer width but after this layer 
impinges on the wall boundary layer, it becomes a constant fraction of 
the duct radius or width. 
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READ INITIAL CONDITIONS 
& INITIALIZE CONSTANTS 


CALL RADIUS I 


SPECIFIES THE 
DISTRIBUTION 
NODES IN THE 
DIRECTION 


CALCULATES THE 
■ DIMENSION LESS DUCT 
'WIDTH OR RADIUS BASED 
ON INPUT DUCT GEOMETRY 


* 


CALL PS I 


CALCULATES THE DIS 
TRIBUTION OF NODES 
Y -PLANE FROM THE * 
DISTRIBUTION 


IN 


CALL PROF 
d CALL YDIS 1 


ASSIGN VELOCITY & 
TEMPERATURE ON EACH 
NODE POINT ACCORDING 
TO INITIAL CONDITIONS 


r- fNSTEP=NSTEP+l| 



CALL LOOkT^I 


CALCULATES JET SHEAR 
LAYER WIDTH & BOUNDARY 
LAYER THICKNESS 


CALCULATES DUCT 
WIDTH OR RADIUS 
AT THIS NEW 
POSITION 


[CALCULATE EDDY VISCOSITY! 


-{ASSIGN STEP SIZE AX I 

, T 

r ~ ■ ■ -fCALL RADIUS! 
(ESTIMATE AP[ 


INNER LOOP 

CALCULATES VEL. K TEMP. ON 
M-2 LINE 


Id 


CHECKS THE 
CONSERVATION OF 
MOMENTUM 


I CALCULATE FRICTION VELOCITY] 

T 

^ | CALL MCHECKl 

- ^ ‘ * 


REPLACE M=1 LINE VALUES 
WITH M=2 LINE VALUES 

I 


LA. 


[call CHECKl^j 
- [print out results] 


CHECKS THE CONSERVATION) 
OF MASS & STAGNATION 
TEMPERATURE 


TERMINATION 
[CONDITION SATISFIED^ 

'yes" 


ZZfj^s 

[stop] 


Figure D-l 

Computer Program Flow Chart 








FLOWCHART - INNER LOO? 



Figure D-l (continued) 





Section (5): Calculation of Velocity on m = 2 Line 


(a) The duct radius or width at the m = 2 line, is interpolated from the input 
data. 

(b) Initially, the m = 2 line pressure gradient is set equal to the average of 
the pressure gradients at the previous two m lines. 

(c) The distribution of velocity on the m = 2 line is calculated. 

Section (6); Calculation of Temperature on m = 2 Line 

(a) The distribution of temperature on the m = 2 line, is calculated, and from 
it the distributions of density and molecular viscosity. 

Section (7): Pressure Gradient Modification 

(a) The position of the nth node, in the y-plane, is calculated from the m = 2 
line profiles. If this value is acceptably close to the duct wall, the pressure 
is incremented by dp. 

(b) Alternatively if this requirement is not satisfied, then an improved 
estimate of the pressure gradient is made. 

(c) Using this estimate, section 5(c) and section (6) are repeated. 

Section (8) : Transference 

(a) The values of u and 9 calculated on the m = 2 line are transferred to the storage 
space previously used for conditions on the m = 1 line, in preparation for the 
advance to the next m-line. 

Section (9): Output 

(a) The velocity and temperature profiles are printed out at defined intervals, 
and several flow variables are printed at every step. 
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Section (10): Termination Test 


(a) If the maximum x-value has not been reached, execution is returned to 
Section (3), in order to advance to the next m-line. The functions of 
each sub-routine are as follows: 


CALC: This evaluates u and 6 using the Thomas algorithm. 

RADIUS: The duct shape is inputed to the calculation procedure, through 

this routine. It interpolates this data and calculates the local 
duct radius at every m-line. 


TEMP: If the dimensionless temperature variation is a known boundary 

condition , it is specified in the routine. This routine is redundant 
with the present boundary conditions. 

YD IS: The position of the grid nodes in the y-plane is calculated with this 

routine . 


PSI: The position of the grid nodes in the ^-plane is assigned in this routine. 

The initial flow conditions determine the form of this routine, i.e, 
single stream flow, two-stream and mass ratio. 

LOOK: The shear layer and boundary layer width are calculated using a 

search technique applied to the m = 1 line velocity profile. 


CHECK: This routine checks the conservation of mass and energy. 

MCHECK: This routine checks the conservation of momentum between adjacent 
m-lines . 


PROF : Calculates the initial velocity and temperature profiles. 
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FORTRAN SYMBOLS 


MEANING 


A(I) 

B(I) 

BB 

BE 

BH 

BY 

CC 

C(I) 

D(I) 

DELTA 

DPI 

DP2 

DP11 

DX 

E(I) 

ENERG 

FDUCT 
FPRIM 
GAMA 
IF LOW 


ITER 

JFLOW 


Minimum value of step size AX 
Dimensionless jet shear layer inner edge 
Dimensionless jet shear layer width, 
BY-BE 

Dimensionless jet shear layer outer 
edge 

K 


c 


n-1 



din 


f at m 
at m 
^at m 


AX 

E 


1 line 

2 line 
0 line 



0uP*AY 


Mass flow in duct 
Mass flow in nozzle 


7 

Control variable (zero upstream of point 
where wall boundary layer and shear layer 
meet, otherwise one) 

Iteration counter for inner loop 
Control variable with the value one for single- 
stream flow and two for two-stream flow 
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FORTRAN SYMBOLS 


MEANING 


LVH 

Dimensional local velocity head, pu / 2g Q 

LZ 

L 

m 

N 

Total number of node points on each m-line 

, NP, NPP, SQP 

Control variables for axisymmetric flow NL 
NP= 2, NPP= 0 and SQP= 0.5; and for plane 
flow NL= 0, NP = 0, NPP=1 and SQP = 1 

NJ 

Number of node points in jet 

NSTEP 

Number of downstream steps 

NTEST 

Test number 

PA MB 

Ambient pressure 

PCUM 

Local dimensionless pressure 

PE 

Pressure at nozzle exit plane 

PH20 

Local dimensional pressure 

PS(I) 

* 

*n 

PR 

P 

r 

PRT 

P rt 

PSN 

$ 

Total in the duct 

PSNJ 

* 

Total ip in the jet 

RHO(m, I) 

* 

P 

RM 

M. 

lr 

RNU 

% 

ROREF 

P o 

RR(I) 

Duct width or diameter 

T(m, I) 

9 

< 1 

TCLI 

T. 

3 


N 

r— 

TFLOW 

Total mass, flow rate, \ p* u AY 


i= 1 


TREF 

TSEC 

TTP(J) 

TTT(J) 


TWREF 

U(m, I) 
UCLI 

UPOT 

URR 

USE 

UUU(I) 

VHEAD 

VIS(I) 

X 

XX 

XX(I) 

XRO 

Y(I) 

YJ 

YS 


Temperature of secondary flow at nozzle exit 

plane 

T 

Dimensional stagnation temperature at each node, 
—2 


T + 


Pu 

2 Sr 


wr 


u 

u 

o 

Velocity at inflection poiiii 
* 

u 

Secondary velocity at nozzle exit 
u 

2 

Reference velocity head, P q u q / 2g Q 

* 

P 

X 

Di stance from duct inlet at which calculation 
ends 

Distance from duct inlet at which duct width, 
RR(I) are provided 

Non-dimensional downstream distance with 
respect to the initial duct half width or radius 
Y 

Half nozzle width or radius at nozzle exit 
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Definition of the Input and Output Parameters 



The input data to the program must be prepared according to the following 

sequence : 


Card No. 


Parameters 


Format 


1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 

14 

15 

16 

17 

18 
19 



Card 19 is the last card to end the calculation of the Program, on which 
NPP must be set a value larger than 1. 

Cards 1 through 18 are required for each set of data. For data more than 
one set, cards 1 through 18 must be repeated in the same sequence. 
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An example of input data for two sets of data are shown on Table A-l. The 
input parameters are: 


SQP , NP , NPP , NL 


DPI, DP2, DP11 


X- 

XX- 

MK- 


PROFE(I)- 


NTEST- 

POl- 

TOl- 

PAMB 

TOO- 

AMASS1- 

AMASSO- 


Control Card for axisymmetric flow 

SQP = 0.5, NP = 2, NPP = 0, NL = 1, for plane flow 
SQP = 1.0, NP = 0, NPP = 1, NL = 0. 

Initial guessed dimensionless pressure gradients on m=l, 2 and 
zero lines respectively. The initial guesses of the values of DPI, 
DP2, DP11 at the initial plane may be assumed equal at any plus 
or minus dimensionless value of the order of 10~ to 10 . 

Distance from duct inlet to nozzle exit plane at which calculation 
begins, inches. 

Distance from duct inlet at which calculation stops 
A number which indicates the number of velocity and 
temperature detail being printed out. 

An array contains MK value of downstream positions in inches 
where the velocity and temperature detail are required to be 
printed out. 

Test or Run number identity 

Stagnation pressure of the primary flow, psia 

Stagnation temperature of the primary flow, °R 

Ambient pressure (i. e. , stagnation pressure of the secondary 

flow), psia 

Stagnation temperature of the secondary flow, °R 
Primary mass flow rate. 

Secondary Mass flow rate. 

For axisymmetric flow, lbm/sec. 

For plane flow, lbm/seo-in. 


RD- Half duct width at nozzle discharge plane, inches 

YJ- Full jet width, inches, (Outside dimension of nozzle exit) 
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Input Data Example for Runs 3 and 6 



PS(I)- To take care of the boundary layer problem, the last 18 values to 

the wall are required to be specified in the SUBROUTINE PSI . 
SUBROUTINE PS(I) already includes the values needed for the 
computer calculation. The computer values were selected to 
satisfy the following: 

(1) grid spacing of the wall should not correspond to a value of 
y + greater than 3. 

(2) neighboring grid spacings should not differ in size by more 
than 50%. 

(3) close spacing is also required in any region away from the 
wall where the velocity gradient is large. 

NS - Indicates the total number of duct geometry to be read in the 

SUBROUTINE RADIUS 

RR(I) - An array contains the total number (NS) of duct width, inches 

XX(I) - An array contains the total number (NS) of axial downstream 

distance, where RR(I) are provided, inches 


Part (b) Output Parameters 


The first section of the output repeats the most important input data for 
different test or run number. 


Velocity ratio 


initial velocity of secondary flow 
initial velocity of primary flow 


Width ratio 


initial duct width 
nozzle width 


Mass flow ratio ~ 


secondary mass flow 
primary mass flow 


J - Indicates node point counting from centerline to wall 

Y(J) - Dimensionless Y coordinate with respect to the local half 

duct width. 


U(J)- 


Dimensional velocity on Jth node, ft/ sec 



TO(J) - 

Dimensional stagnation temperature on Jth node, °F 

I - 

Print step counter at approximately XIN increases 
0.25 inches 

XIN - 

downstream distance from nozzle exit plane where 
calculation begins, inches 

X/BO - 

ratio of downstream distance with respect to initial 
half duct width. 

B/BO - 

ratio of the local duct width with respect to initial duct 
width. 

PH20 - 

Local wall static pressure, 9 inches of water 

UCENT - 

velocity of the flow at centerline, ft/ sec 

TOCENT - 

centerline stagnation temperature of the flow , F 

y 

AUGMENT - 

J W 

> w _ _2 

Local momentum flux, 2 / p u dy , divided 


o 

by initial jet momentum 


USTER - Local friction velocity, ft/ sec 

The selection of intervals at which calculations are made is determined by 
a subroutine in the computer program. The data is printed out at approximately quarter 
inch intervals. The locations where temperature and velocity profiles are printed out 
are specified by the user in PROFE(I) described in the input data. 
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PROGRAM LISTING 


The program listing included in this report is for the program as run on a 

CDC/6600. The program was initially developed on an IBM 360/50. The deck is the 
one successfully run on the CDC/6600. 

The essential changes to the program necessary to recover the IBM 360/50 

deck are: 

I. Certain variables should be in a real* 8 mode 

Add 


REAL* 8 Y, DABS, DLOG, YB1, YB2 
REAL* 8 ZY, Y 
REAL* 8 Y 
REAL* 8 Y 
REAL* 8 Y 

Functions ABS( ), and ALOG( ), should be changed to DABS( ), 
and DLOG( ). These occur in cards 

JMX2710 JMX3750 JMX4450 

JMX3030 JMX4340 JMX4460 

Card YIS50 should be Y(l) = 0. 0D0 

Output Hollerith symbols should be changed from * to These occur in 

cards 


JMX460 

JMX660 

JMX730 

JMX2120 

JMX5640 

JMX470 

JMX700 

JMX740 

JMX5190 

JMX5650 

JMX520 

JMX710 

JMX2090 

JMX5200 

JMX5660 

JMX530 

JMX720 

JMX2110 

JMX5590 

JMX5670 

GEM110 

CHE330 

BLC190 



GEM120 

CHE340 

BLC200 




Deck characters are BCD. 
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Figure D- 2 

COMPUTER PROGRAM LISTING 


RUN VfPSTON Z.7 — PSR LEVEL 29*-- 


OOOOOI 

000003 

OOOOO? 


000003 


c 

c 

c 

c 

c 

c 

c 

c 

c 

r 

c 

c 

c 

c 

c 

c 

r 

c 

c 

c 

c 

c 

c 

c 

c 

c 


ooool* 
0 00 1 6 
00001 * 
oocnl* 
00003? 
00003? 
000035 
000041 

000041 
000043 
0 0 0044 
oooo4<; 
0000*51 

oooas i 

OOOOS-5 

000054 


PROGRAM NAS (INPUT, OUTPUT, T APE5* I NPUT ,T APE6-0UTPUT) 

PEAL LZ 

COMMON NP ♦ 5GP * NPP t NL % 

DIMENSION 11(2,70 l.miJUI. T(2«70 ) tSl (70 >,S?(70 ),S3(70 )• 

1 PS ( 70 ) , R2 ( 70 > .El (70 I.PK70 ), 7(70 ).E{70 )» 

3PHO(?*70) .VIS(7a> ,TTP(70) *AA(70> » 0(701 
T,A(70 ).R(70 ) »C (70 )»H(70 )»BBB(70) 

4 ,Tn L 0(2> ,AUGM(70) ,PR0FE(20) *UJ(70) .7RL(701 
F ( X 1 *Y1 ) 3 ( ( 2 , /? • A ) * (1 ) ** 3 “7(*X1 

• * • SECTION i • * • 

K 1 TS A CONTROL VARIABLE. 1 EOR SINGLE STREAM FLOW WITH I^ITI 
rOI iNpARY LA YER . 2 FO* SINGLE STREAM FLOW WITHOUT W.B.L. OR TW 
FLOW WITH A TOP HAT PROFILE 

X *0lSTANCE FPnM THE DUcT inlet to nozzle exit INCHES 

^xx st ance^fpom ^ the 5 duct INLET AT WHICH CALCULATIONS stop 

PFs*OZZl| R |mt PLANE S HRES$URE LBF/FT2 GAUGE 
ROPFF=REFERENCF DENSITY LPM/FT3 
FOUCTaTOTAL MASS FLOW RATE ILBM/SEC ) 
ppspHANTL NUMRfR 
PPT = TURBULENT PRANTL NUMRE& 

T wrFF*W ALL REFERENCE TEMPFrATURE OEGr*R 
TCLI= JET TEMPERATURE AT NOZZLE EXIT OEG*R 
YJs EFFECTIVE nozzle EXIT RADIUS t INCHES) 

UCLT * JET VELOCITY AT NOZZLE EXIT FT/SEC 
GAMAsGAS CONSTANT^ 1*4 FOR AIR 
RMa REFERENCE MACH NUMRER 
TSFC = SECONDARY TEMPERATURE (DEGr*R) 

PRFF=REFF«ENCE Pressure LRF/FT2 A 
PAMRsAMRIFNT PRESSURE LRF/FT? 

N07ZLE RADIUS Im INCHES 

OAT A RNU, ROREF t PREF, GAMA/, nOOl 590 * . 0763 , 21 15 . 1 1 *4/ 

OAT A PR, PR Tt TREE* TWREF/oWt 0*9 t520«0f 560*0/ 

RA89 CONTINUE 

PFAD(5»8olR> SOPtNP i NPP*NL 
p^Iq F OR M A T ( F 5 t 0 ♦ 3 1 5 } 

IF(NPP- 1) 737,838,883 

?79 format mhi»//i iox»*two-dimensional ouct 

*.//) 

TWRtF »54 0 • 

CC ■ .08 
GO TO 1023 

7 fl2q FOPMAt(lHU//.lnX, .AXISYMMETBIC DUC7 FLOW* , / , l OX .*— -— ** 

t ».//> 

TWRfF a 560. 

CC = .07 

r INPUT SECTION 

CONTINUE 

1 


JMX 

50 

JMX 

60 

JMX 

70 

JMX 

80 

JMX 

90 

JMX 

100 

JMX 

110 

JMX 

120 

JMX 

130 

JMX 

140 

JMX 

150 

IAJMX 

160 

’R JMX 

170 

JMX 

180 

JMX 

190 

JMX 

200 

JMX 

210 

:hjmx 

220 

JMX 

230 

JMX 

240 

JMX 

250 

JMX 

260 

JMX 

270 

JMX 

280 

JMX 

290 

JMX 

300 

JMX 

310 

JMX 

320 

JMX 

330 

JMX 

340 

JMX 

350 

JMX 

360 

JMX 

370 

JMX 

380 

JMX 

390 

JMX 

400 

JMX 

410 

JMX 

420 

JMX 

430 

JMX 

440 

JMX 

450 

.--JMX 

460 

JMX 

470 

JMX 

480 

JMX 

490 

JMX 

500 

JMX 

510 

-- JMX 

520 

JMX 

530 

JMX 

540 

JMX 

550 

JMX 

560 

JMX 

570 
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R, t N v^PSTON 2.1 — P SR LFvFL ?9 ft-- nA S 

00005* READ ( 5 * 7R 1 7 ) DP1,DP2,DPU 

000066 7R17 FORMAT ( 3E13.6) 

000066 RF An ( 5 * 73) x,XX 

000076 RF *0 ( 5 » 17 ) MK 

RFA0(5f 73) CPROFE(I) *t*l ,mk> 

0001 1 7 RF AD ( 5 1 1 7 ) NTFST 

UOUl?* 37 F0RMATU5) 

000125 WPITE(6,30) NTFST 

000131 19 FORMAT r//,25X f ; INPUT DATA RUN NO. ♦,l5 t //> 

5°?! 11 , HF AO (5.73) POl ,T01 ,PAmP # TOO* AMASSl , AMASSO, RDtYJ 

000157 73 FORMAT (RF1 O.nj 

0JJO157 WRJTE (6.7) POl ,T01 .PAMR.TOO, AMASSl ,AMASSO,ROtYJ 

000203 A 7 FORMAT ( 1 5X » *PO 1 a *.Ffl.4,* PS I A* f / , ] Sx * *TO 1 • *»Ffl.4 # * 
l/t15Xt*PAMM a # * Ffl , 4 * » PSIA*,/t 15Xt«TOO a *,FB. 4 t* DEG, 
?/* l5Xt*AMASSl a *,F8,4,* LRM/S-IN ,* * / f 15 X, »aMASS 0 * 
3/S-IM.*t/tl5Xf#R0 (HALF DUCT WIDTH) a **FB,4 1 * IN.*./»15X. 

4 JFT WIDTH) o * f Ffl.4,« lN.*./> 13 ' 

000PQ3 I F l NP .fO, 2 ) r .O TO 20M 

00020R A? * 2 * *R0 -YJ 

000210 GO TO 1201 

000211 2001 A? * 3, 14) 6/4, * { (2 # *R0) **2-y J#*2) 

0O0216 1201 CONTINUE 

0t ' 0?1 * ASTApa{ANAssO/PAMB)*$op T ( T 00* < 5 3 . 3/ ( 1 ♦ 4#32 . 2 ) ) )*1.2**3 

000224 Y 1 x A2/ASTAR * 

000230 Xll s 1. 

00023) X?2 » ,03 

000233 TDt a ,0001 

000 234 TE ST 1 *F ( X 1 1 t Y 1 ) 

000237 TEST?*F (X22» VI) 

000241 ?? HALF 8 ( X U ♦ X 22 ) /2, 

000244 CHFK1 =F {HALF 1 Y1 ) 

>00247 IF (AF^S (CHEK1 ) .LF,TOL) GO TO 21 

300252 TF { ARS { TEST l-CHFKi)-ABS( TEST T) j 19,23*20 

)00?56 ?9 X2? a HALF 

>00260 TEST2=CHEK1 

>00261 GO TO 22 

>00262 19 XI 1 * HALF 

)0 02 64 TEST) eCHEK) 

10 0 264 GO TO 22 

100266 21 RM? = HALF 

,0 0?7f> PI a Pa^B/ ( ( 1 , ♦ ,p* RM ? «*->) 

'00276 PA1 e (P01/P1 ) «« ( { G AM A- ) ) /G AM A ) 

>00305 Rm = SORT { ( P A 1 •],)*(?«/( GAMA«i») ) ) 

'00 31 A TcLT * T01/PA1 

'00116 UCL I a RM*SURT ( 1 ,4»32,2*53,3«TCUI) 

'00321 P A2a { P A W R /P 1) »# ( (GAMA- ] . ) /GAMA) 

'00132 TSEC * TO0/PA2 

'00334 gSFC * SORT f?,*.24*778.*32.2* (TOO- TSEC) ) 

00141 PF e (Pl-PAMB) *144, 

00344 PaWrs P a Mr * 1 4 4 , 

00346 U JR a USFC/UCLI 

00347 AMR a AMAS50/AMASS1 

00351 YJ a YJ/2, 

0035? KJa70 

00351 NTT=N-2 


JMX 

5*0 

JMX 

590 

JMX 

600 

JMX 

610 

JMX 

620 

JMX 

630 

JMX640 

JMX 

650 

JMX 

660 

JMX 

670 

JMX 

680 

JMX 

690 

♦ JMX 

700 

JMX 

710 

iM JMX 

720 

L JMX 

730 

JMX 

740 

JMX 

750 

JMX 

760 

JMX 

770 

JMX 

780 

JMX 

790 

JMX 

boo 

JMX 

BIO 

JMX 

820 

JMX 

840 

JMX 

850 

JMX 

860 

JMX 

B70 

JMX 

880 

JMX 

890 

JMX 

900 

JMX 

910 

JMX 

920 

JMX 

930 

JMX 

940 

JMX 

950 

JMX 

960 

JMX 

970 

JMX 

9B0 

JMX 

990 

JMX 

1000 

JMX 

1010 

JMX 

1020 

JMX 

1030 
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RUM VT0SIOM ?.3 — PSR level 298— 


nAS 


00*35* 

0003*7 

000^60 

000161 

00016? 

000361 

000364 

00036* 

000366 

000167 

00037P 

000371 

00037? 

00037? 

000174 

000175 

000176 

000377 

000400 

000401 
000404 
000406 
000407 
00041 ? 
00041 * 
000417 
0004?0 
0004? 1 
0 00 4 ? 1 
0004?4 
0 0 0 4 ? 6 
00043? 
000437 
000444 


000444 

00045* 

00045? 

000454 

000455 

000460 

000463 

000466 

000471 

000474 

000476 

000*00 
000501 
OOO5O7 
0 0 05 1 0 


C 


1* 


NNa N-l 
PPFF « 2 1 1 5« 
PT =• 0. 

TD| 0(1) « 0 • 
ITFP a 0 
BVsn, 

IL* 0 
JFLOW«2 
UPFF * t 00l 
IURF.S»0 

BHsO, 

ph?o * n. 

N S T F P B 0 
IFLOW a 0 
URB * l. 
TCOrEsO 
TOLIH?) *0. 
DFL TA = 0 • 


KJ*? 

UJBsU SECOnDARY/U primary 
PEN*UCL 1 / ( ?4 • 4RMIJ) 

IF <NP .EO. 2) GO TO 118 
AmASSI » AmA$$1«12. 

FOUCT a AMA 550* 1 2 • ♦ AMASSl 
PSN = FOUCT/(?.*RNU*flORFF) 
pSNJ=AMASS) /(2.*RNU*R0RFF) 


FOt'CT * PSM 
GO TO 117 
CONTINUE 

YOUT 3 RD*?**PFN 
FOuFy = AMASSl ♦AMASSO 
PSN sFOUCT^UCLI/ (6.284*RNU*42*P0R£F) 
pcKj 1 K AmASS1*uClI/ (6,?84*RNU**2*RQREF) 
FDuCTaEOUCT*OCL T / { ROPEF# 3 . 1 4 2#YOUT**2*HNU # *2 ) 


117 CONTINUE 

SPECIFY STREAM FUNCTION DISTRIBUTION 

call PSI< n, psn,nj,ps,psnj> 

00 R J=?*NN 
JP 1 * J ♦ 1 
J M 1 a J — l 

51 (J)=PS(JP) )-P5(J) 

52 ( J) =P5( J) -PSf JM1) 

53 (J) =PS <J p l ) “PS (JM1 ) 
b 1 (J)=S1 (J)/S2(J)/S3(J) 
B2(J)aS2(J)/Sl(J)/S3<J> 

R CONTINUE 

no 9 I a 1 t N 

<5 FIT) a 0* 

X *x *UCL I / <RNU*1 2 • ) 
XXsXX*UCLI/<RNU*12,) 
yJ*YJ4uCLI/(HNu*12.) 


JMX 1130 
JMX 1140 
JMX 1150 
JMX 1160 
JMX 1170 
JMX 1180 
JMX 1190 
JMX 1200 
JMX 1210 
JMX 12?0 
JMX 1230 
JMX 1240 
JMX 1250 
JMX 1260 
JMX 1270 
JMX 1280 
JMX 1290 
JMX 1300 
JMX 1310 
JMX 1320 
JMX 1330 
JMX 1340 
JMX 1350 
JMX 1360 
JMX 1370 
JMX 1380 
JMX 1390 
JMX 1400 
JMX 1410 
JMX U20 
JM X 1425 
JMX 1430 
JMX 1440 
JMX 1450 
JMX 1460 
JMX 1470 
JMX 1480 
JMX 1490 
JMX 1500 
JMX 1510 
JMX 1520 
JMX 1530 
JMX 1540 
JMX 1550 
JMX 1560 
JMX 1570 
JMX 1500 
JMX 1590 
JMX 1600 
JMX 1610 
JMX 1630 
JMX 1620 
JMX 1640 
JMX 1650 
JMX 1660 



RUN VFRSTON 2.3 --PSR LEVFL ?9fl 


NA5 


00051? 
000s 1 6 
000517 
0005?! 

0005?! 

000524 

000526 

000531 

000535 

000540 

00054.3 

000547 

000550 

000551 

000557 

000555 


000556 

000566 

000571 

00057? 

0005 7 4 

000576 

000507 

000610 

000621 

000423 

000641 

000645 

000647 

000653 

000653 

000666 


000666 


000667 


000670 


CALL RAOIU$(X,YOUT, 1 ,REN> 

YTsYoUT 
NTP=N-3 
81 FNO a 4 ,*YT 

c hi.fnuing length in radii 
xBLFND=0 , 

XHO=X/YOUT 

c ♦N*TI>>LI5 aTI0N of CONSTANTS 

IF OiJk.G7.,990)JFL0w« 1 
IF f JFLOrt »EO« 1 )NJ*N 
IF f JFLOW.fTO. 1 ) cr*0.09 
C INITIALISATION OF QP/DX VARIABLES 

VHFAO = UCL I4* 2 *^riPEF/64*4 
T (TSFC-yCLI ) /( twrEF-tCLI) 

PHFFsPWEF /VHEAD 
PAviRsPAMR /VHEAD 
PCtiMaPE/VHEAD 
NNaN-l 
URP * 1, 

C 

c 

c 

C * * * SECTION 2 * # # 

C 

c *♦******##*•**## 

C INITIAL FLO* CHARACTERISTICS AT Mxj LINE. 

C VFLOCITY AND TEMPERATURE ON Mai LINE 

CALL PROF <U» T *Nj,N,UJ8,TJR,KJ,UCLI) 

C 

CONAs (TWREF-TCLI ) /tref 
COnP = TCl T /TREp 
CONFxl9B,7/TPEF 
rONpal .♦CONE 

CONI. 3( r,AMA.i )# RM **? ./(TyREF/TCtUl.1 

no 5 t 1*1 tN 

RhO n * I ) = (P C UM*P ft M B ) /(PP FF « ( C 0N A *TU t J> ♦rONftl ) 

RHO(? f l)a«HO(l f Il 

61 V 1 6 { i ) = ( T ( i tI)#C0NA*C0Nfn##1.5»C0ND/(T(lf U*CONA4CON8*CONE) 
r 

CALL YDI$(YtP5,RHOfU»N) 

OJP=YOUT/YJ 
WPtTE<6, lo51) 

1051 FOPMATdHl,* .PRINT OUT INITIAL VALUES*,//) 

WRITE (6, 7232) I.IJR,OJR,Amr 

7232 FPRMAT(20X,«VELOCITY RATIO * * i F 1 4 . 4 f / * 20 X , *W I DtH ( 0 1 A) RATIO 
1F14.4,/,20X ,*MASS FLOW RATIO * *.F14.4,/) 

C JJK IS THE NUMBER OF THE GRID NODE AT YJ PLUS 2 
' JJK*NJ*2 

C * * * * * * * 

C PRINT OUT INITIAL VALUES 

C 

GO TO 1748 

C ******* 

C INITIALISATION E^OS ♦OUTER 

loo CONTINUE 


* * * 

LOOP BEGINS 


JMX 1670 
JMX 1680 
JMX 1690 
JMX 1700 
JMX 1 7 1 0 
JM X 1720 
JMX 1730 
JMX I74O 
JMX 1750 
JMX 1760 
JMX 1770 
JMX 1780 
JMX i 79 0 
JMX 1800 
JMX 1810 
JMX 1820 
JMX 1 830 
JMX 1840 
JMX 1850 
JMX 1860 
JMX 1870 
JMX 1880 
JMX 1090 
JMX 1900 

* JMX 1910 
JMX 1920 
JMX 1930 
JMX 1940 
JMX 1950 
JMX I960 
JMX 197 
JMX 1980 
JMX I99O 
JM X 2000 
JMX 2010 
JMX 2020 
JM X 2030 
JMX 2040 
JMX 2050 
JMX 2060 
JMX 2070 
JMX 2080 
JMX 2090 
JMX 2100 

a * » JMX 2110 
JMx 2120 
JMX 2130 
JMX 2140 

* JMx 2150 
JMX 2 160 
JMX 2170 
JMX 2180 

* JMX 2190 
JMX 2200 
JMX 2210 
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RUN VERSION 2.3 — P SR LEVEL 29A — 


NAS 


000670 

00067? 

000671 

000674 

0007U 


000711 
00071 1 

00071* 

0007)7 

00072O 

00072? 

000721 

0007?* 

000726 

00073* 

000737 

00074? 

000744 

000755 

000757 

000760 

000760 

000764 

000777 

00100 ‘s 
001 00 6 
001007 
00101) 

001011 

ooion 


001044 

001047 


NSTEP*NSTEPM 

c 

c calculate fhictton velocity 
c 

lF(NSTEP.EO.i)GO TO 631 

c 

c 

IF ( IFLOW ,E0* 1 ) GO TO 80 

c 

CALL LOOK( JJK,NN,UtY*DELTA f 6 H , 0Y * N # CC * IFLOW * XBLEnD* Y J» BE 1 

c evaluates smear layer width ecT f 

C 

*0 CONTINUE 

c SECTION 3 INSERTED HERE 

C * * * SECTION 3 * # # 

c 

C * EVALUATE TURBULENT VISCOSITY USING NIXING LENGTH # 
c VERSION for pipe flow (AS PER APRIL 17* 1972 ) 

C IF ( I FLO* , £Q , 1 ) GO TO 770 

tEUFLOw.EO.nOO TO 77ft 
C jFLOw * 1 INDICATES A PIPE FLOW 
*DELT a »Y(N)-OELTA 


1*2 


TNM*0 

A2l ■ *oR*OEXTA 
DO 4* 1*2* NN 
TP1 *t ♦ 1 
I Ml *| -1 

DELUaAHS ((II ( 1 ♦ IPO -u <1 * TM1 )) /U ( 1 * I > > 

IF(0FLU.LE.. OOlOGO TO 1360 
IF(Y(I).LT.0E)GO TO 1360 
IF ( I *EO# 2) GO TO 796 

IF POINT IS BETWEEN SHEAR LAYER AND WALL L » GO TO 1360 
IF (Y ( IMl) .GT.GY.ANO.Y (I) .LT.ROE lT AjGO TO 1360 
IF POINT IS IN WALL 0. L. GO TO 1 3 T3 
IF(Y(I) .GE.BOELTAJGO TO 1373 
GO TO 796 
WALL viscosity 
CONTINUE 

YS*(Y(N)-Y(I) )*UREF 
L Z*0.41*(1 # 0*ExP(-YS/26.0))*IY<N >-Y<I >> 

IF ( L7 • L T • A71 ) INm*INMM 
IF (LZtGT. AZDLZsAZl 
GO TO 76? 

1360 E d ) *0 * 

GO TO 45 
JET VISCOSITY 

F(I > * A0SUY II) -Y (IMl ) )*<Un •IP1> *UM * I) 1 / MY (IP) ) -Y ( H )• 

’ (vnpn-YdMn ) )♦ (Y(ipn-Y(tn*(un *n-un *iwn >/( (Ydj.YdMi) 

?<Y( IPI)-Y(IMI) n ) *L2*LZ 


1173 


796 

76? 


45 


CONTINUE 
GO TO 48 
PIPE FLOW 


MIXING length 


JMX 

JMX 

JMX 

JMX 

JMX 

JMX 

JMX 

JMX 

JMX 

JMX 

JMX 

JMX 

JMX 

JMX 

JMX 

JMX 

JMX 

JMX 

JMX 

JMX 


JMX 

JMX 

JMX 

JMX 

JMX 

JMX 

JMX 

JMX 

JMX 

JMX 

JMX 

JMX 

JMX 

JMX 

JMX 

JMX 

JMX 

JMX 

JMX 

JMX 

JMX 

JMX 

JMX 

JMX 

JMX 

JMX 

JMX 

JMX 

JMX 

JMX 

JMX 

JMX 

JMX 


2220 
2230 
2240 
2250 
2260 
2270 
2200 
2290 
2300 
2310 
2320 
2330 
2340 
2350 
2360 
2370 
2300 
2390 
2400 
2410 
2420 
2430 
2440 
2450 
2460 
2470 
2400 
2490 
2500 
2510 
2520 
2530 
2540 
2550 
2560 
2570 
2500 
2590 
2600 
2610 
2620 
2630 
2640 
2650 
2660 
2670 
2660 
2690 
2700 
2710 
2 720 
2730 
2740 
2750 
2760 



vroslO*' 2.7 — PSR LEVEL 298 — 


NAS 


-jn l *4 7 
on l "50 
nn i n51 
00 1 "53 
on 1 ns* 
90)057 
oninAi 
OO 1 n6? 
001*67 
00 1 *66 
001073 
001076 
001076 
00110 ] 
oono* 
001107 
001116 
001116 
00113? 

oom* 

ooiuo 

001153 

00115* 

001160 

001161 

001163 

001165 


00i?16 

001221 


001221 


001221 

001223 

001226 

001231 

001234 

001236 

0012*2 

0012*6 

001252 

00)256 

001262 

001263 

001265 

001267 


770 TNM.-i 
IT*0 

AZ * .09.0ELTA 
00 In J*2*NN 

IP1 * I ♦ 1 
IM1*I-1 

IF ( TFLOW ,E0 • 1 ) GO TO *172 
UNN«,00604|J(1, 1) 

UN*AftS(UU*THl)-UU*IPlM 
IFCUN.GT.UNN > r,0 TO 81*0 
GO TO 10 

**1 72 Yn.Y(N)-YU) 

yd=i t-yn/v{N> 

TF<X8LEN0.LT.BLFN0> GO TO *173 
AZ* ( . U-*^e*Y04»2-.06*Y04#4) *Y (N) 

GO TO B 1 40 

4173 AZ *CC *PM*{ (.l4-*084Yn*«2- #0 64YO*44)*Y<N)-CC # flH) *XBLENO/0LEND 
81*0 TFnt.EO.DGO TO 832 
YS«(Y(N)-YU))*UREF 
L7 ».4l*( U-EXPI-YS/26.) ) # (Y(N)-Y(in 
I NM* I NM ♦ l 

IF(L7#GE,AZ) IT*] 

AZ1«lZ 

TFCLZ.LT, A7) GO TO \\ 
fl3? AZ1*AZ 

11 F( T)* ARS( (Y(I).Y(TM1))*(UC1*IP1)-U(1,I) )/{ CY(IPl)-YcI))# 

1 ( Y ( IP1 ) -Y ( IM 1 ) DMY(IPn-Y(I) ) # <U(l*I)-U{lfIMD )/< ( Y ( I ) -Y { IM1 ) ) # 
2 ( Y t I p 1 ) -Y ClMl ) ) ) ) *AZ 1* AZ l 
10 CONTINUE 
*8 CONTINUE 
C 

C FINAL CARD OF SECTION 3 

631 CONTINUE 

C SECTION A INSERTED HERE 

C # # # SECTION * * # • 

C 

C EVALUATE Ox STFP <FOR EJECTOR FlO# AS PER FEB# 19) 

OX»Rh* ,* 

BB».0*04YCN) 

IF(DX.LT.RB)OX*PR 
7K*YCN)/8. 
lF<r>X.GT.ZK)0X*7K 
fF (NSTEP.LT. 10) nx*,54QX 
IF (NSTEP.LT #20) 0X*,2*DX 
IF (NSTEP.LT. 30) DX«,5*0X 
TF ( IuRES.Eo, 1 ) 0)(*,3*0x 
TF(XWO,GT.]2,0) OX*0X*2. 

IUBfS*0 
12 X«X#DX 

XHLEnO*XBLEND4DX 

XR0*X/YT 


FINAL CARO OF SECTION * 


JMX 2770 
JMX 2780 
JMX 2790 
JMX 2800 
JMX 2810 
JMX 2820 
JMX 2830 
JMX 28*0 
JMX 2850 
JMX 2860 
JMX 2870 
JMX 2880 
JMX 2890 
JMX 2900 
JMX 2910 
JMX 2920 
JMX 2930 
JMX 29*0 
JMX 2950 
JMX 2960 
JMX 2970 
JMX 2980 
JMX 2990 
JMX 3000 
JMX 3010 
JMX 3020 
JMX 3030 
JMX 30*0 
JMX 3050 
JMX 3060 
JMX 3070 
JMX 3080 
JMX 3090 
JMX 3100 
JMX 3110 
JMX 3120 
JMX 3130 
JMX 31*0 
JMX 3150 
JMX 3160 
JMX 3170 
JMX 3 lgO 
JMX 3190 
JMX 3200 
JMX 3210 
JMX 3220 
JMX 3230 
JM X 32*0 
JMX 3250 
JMX 3260 
JMX 3270 
JMX 3280 
JMX 3290 
JMX 3300 
JMX 3310 


74 



quN vrOSTHM 2*1 — PS« LEVEL 29A-- NAS 



C 

C * 

• * SECTION 5##* 


C 

C 

determine HOUNDxRY conditions on m#i LINE 

001 ?7A 


U(?tN) *0*0 

001*7? 


vin-Oi 

001271 


CALL RAOIUS(X f YOUTf2*REN) 

0f>l?76 


r)p?«(DPn*0Pi>/?. 

001 3<>1 


I chx *0 

00130? 


TTEP-1 


c 

c 

CALCULATE VELOCITY U ON m*i LINE 

00130? 

t 

■K 

ei m • o. 

001304 


00 7347 J *2,N 

001305 

7347 

El f J) *Y < J ) * # NP *U(1 * J)*PHO(ltJ)/(2**PS(J) 1 

001326 


NTk*NTT-1 


C 

c 

ITERATION STARTS 

001330 

c 

6001 

CONTINUE 

001330 


DO 13 J«2,NTT 

00133? 


JPt*J*l 

001334 


JM1*J-1 

001335 


Y1 *F 1 <JPl)*<VI$tjPl)*RHCHltjPl>*E(JPin 

001342 


Y2«E1 < JMl)#(Vls< JM1) ♦RHO(l,JMlJ#E(JMin 

001347 


Y3*E) ( J) * (VIS { J) ♦RHO (1 tJ) ♦£( J) > 

001354 


Y4* ( Yl ♦ V 3 ) /SI < J) 

001357 


Y5«(Y3^Y?)/S2(J) 

00136? 


Y6*Y4/($3(J>*2.) 

001365 


Y7*Y5/(S3(J) # 2.) 

001367 


Y8*Y6 # U (1 t j)/RS«j) 

00137? 


Y9.Y7*U ( l l J)/PS(Jj 

001375 


A UM1) «U < 1 * J) /0X*Y8*Y9 

00140? 


RUM1)»-Y8 

001404 


C(JM1)«-Y9 

001405 

13 

0 |JMn*U(l»J) *»?/0x-*DP2 / (2,*RHO(l »JJ 1 

001417 


AdltAUXU) 

0014?! 


TF (NP *EQ • 2) On TO 119 

001423 


DIO* < Y (NN) -Y (Nt) / (Y (NTT1-Y (Nil 

00143? 


OBAR«OP1*(Y(nN>** 2-Y <M **?- ( Y<NTT)**2-Y<N)**2)*OlG>/4 # /Vt$<NN) 

001446 


GO TO 219 

001447 

119 

CONTINUE 

001447 


n I <WLOG<Y<NN1 /v <N) ) /ALOG < Y <NTT) /Y CN) ) 

001463 


DBAR«DP)*(Y <NN)**2-Y (N) **?- (Y (NTT 1 **2-Y <N) **2> #0!<j> /B f /V 1S<NN> 

001477 

?1 9 

CONTINUE 

001477 


A(NTK)*A <NTK) ♦r>TG«B(NTK) 

001503 


n<NTK) *D (NTK ) -DRA«*R (NTKl 

001506 


CALL CALC<A*H,C»0,H*NN) 

001513 


00 1934 J*?tNTT 

001514 


U ( ? * J) *H <J-1) 

00I5U 


IF<U(2tJ) .LE..001I IURE3*1 

0015?? 


TFCIJCZt J1 «L£»»00l)U <2fJ)*»00l 

0015?7 

1934 

CONTINUE 

00163? 


ll(?fl>«U(2'2> 


JMX 3320 
JMX 3J30 
JMX 3340 
JMx 3350 
JMX 3360 
JMX 3370 
JMX 3380 
JMX 3390 
JMX 3400 
JMX 3410 
JMX 3420 
JMX 3430 
JMX 3440 
JMX 3450 
JMX 3460 
JMX 3470 
JMX 3480 
JMX 3490 
JMX 3500 
JMX 3510 
JMX 3520 
JMX 3530 
JMX 3540 
JMX 3550 
JMX 3560 
JMx 3570 
JMX 3580 
JMX 3590 
JMX 3600 
JMX 3610 
JMX 3620 
JMX 3630 
JMX 3640 
JMX 3650 
JMX 3660 
JMX 3670 
JMX 3680 
JMX 3690 
JMX 3700 
JMX 3710 
JMX 3/20 
JMX 3730 
JMX 3740 
JMX 3750 
JMx 3760 
JMX 3770 
JMX 3780 
JMx 3790 
JMX 3800 
JMX 3810 
JMX 3820 
JMX 3 h 30 
JMx 3840 
JMX 3850 
JMX 3860 


Bi; m vroSION 2,1 — PStt LEVEL 29S— NAS 


001*3? 


0 (2t NN) »U(??NTT) # OlG*DBAO 

JMX 3870 


C * 

• * SECTION 6 * • * 

JMX 3880 


r 

EVALUATE TEMPFPATUREfVTSCOSlTYf AND DENSITY FIELDS 

JMX 3890 


c 

This section should follow section ? foh incompressible flow 

JMX 3900 


c 


JMX 3910 

00 1 ** n 


DO 1 035 J*?fNN 

JMX 3920 

001541 


JP1 *J^1 

JMX 3930 

001543 


JM 1 a 1 

JM X 3940 

001544 


Yl»ei(JPi>*(VlStjPl>/PR*E«JPl>»ftHO(l*JPll/PRT> 

JMX 3950 

001553 


Y2*El ( JM 1 ) * ( V I S ( JM l)/PP + f( JM 1 ) * RhO (I t JM 1 ) /PRY ) 

JMX 3960 

00 1 5*? 


y3*E \ ( J>* <vIS(J)/PR* (E ( J)4RH0(1 • J) )/P«T) 

JMx 3970 

001570 


Y4» ( Y 1 ♦YB ) /S 1 (J) 

JMX 3980 

001573 


Y5* ( Y3* Y? ) /S2( J) 

JMX 3990 

00157* 


Y 6 *Y*/<$ 1 < J)*2.) 

JMX 4000 

001*01 


y7bY5/(S3(J)*2,) 

JmX 4010 

001*03 


Y8*Y**U<1»J)/P$(J> 

JMX 4020 

00 l*O* 


Y9»Y7*U(1* J)/PS(J) 

JMX 4030 

001*11 


A(jMn»U(l*JJ/DX 4 Y 8 *Y 9 

JMX 4040 

001 * 1 * 


R ( JMi ) *-Y R 

JMX 4050 

001*20 


C<JMl)*-Y9 

JMX 4060 

001*21 

]93* 

n(JMl)3TUtJ)*UUtJl/DX*CONL *0P2 * U ( 1 » J) / < 2« + RHO ( 1 » J) ) 

JMX 4070 



1 ♦ (VTSU1 *E ( J) *«H0(1 * J) ) *PH0(lfJ) *C0NL* 

JMX 4080 



? (M? ( J) # ( U < ? • JP 1 ) -U (2 * J) 1 ♦Ml ( J)*(U(?#J1-U(2* JMl) 

JMX 4090 



3Y(JM*NP *U(lt J)2U(l«Jl/(4,«PSU)®P$(J) 1 

JMX 4100 

001*72 


OIG*l • 

JMX 4110 

001*74 


ROUND* { Y (N- 1 ). YIN -211 4 IV <N) #g. *Y(N-n-Y(N-21)/4* 

JMX 4120 

001704 


0RAqsC0NL*PH4RnUND* <U«EF*V IS I N-l ) 1 »*4*/ (HHO U *N.l > #* 2 .) 

JMX 4130 

001720 


TF(NSTEP*EO,nORAR« 0 . 

JMX 4140 

001723 


A (NTK ) *A { NTK ) ♦0TG«tj (NTK1 

JMX 4150 

001727 


D(NTK) *0 (NJK) •r)RAR*B (NTK) 

JMX 4160 

00173? 


A<l)mA(l)*C(l) 

JMX 4170 

001734 


CALL CALC I A trj »C *D * H *NN) 

JMX 4 I 8 O 

001740 


DO J* J»2|Nj^ 

JM X 4190 

00174? 

1* 

T(?*J)*H < J- X > 

JMX 4200 

001750 


T (? *NN) *T t ? * N t T 1 •OIC^OBAp 

JMx 4210 

00175* 


T(?*N)*T(2*N-l) 

JMX 4220 

00 1 7*1 


T<?*11*T(2*2) 

JMX 4230 

0017*? 


00 50 1*1 * N 

JMX 4240 

0017*3 


RHO ( I * I ) *RH0 (2 • T ) 

JMX 4250 

001767 


RHO (?»!)* (PCUM + pAMB) / (PRf F* ( C0NA*T (2 » I J ♦CONR> ) 

JMX 4260 

00177* 

*0 

Vi S ( I ) * (T (1 .1) •CONA^CONR) ••1,5*C0nD/<TO * I)*C0NA*C0NB4C0NE1 

JMX 4270 


c * 

* * SECTION 7*4* 

JMX 4280 


c 


JMX 4290 


c 

# 4 # * CALCULATE pressure ghaoient 

JMX 4300 

002014 


CALL YOIS (Y*PStPHOiU»N j 

JMX 4310 

002020 


IF(ItE«.GT.9)GOTO 6732 

JMX 4320 

002024 


ITER*tTER*l 

JMX 4330 

0020?* 


YTFST* APS ( YOUT -Y (N) ) 

JMX 4340 

002030 


YPPT*YOUT4*OOOOl 

JMX 4350 

00203? 


TF(YTEST,LT.YPPT >GO TO 6732 

JMX 4360 


c 


JMX 4370 

002034 


tFdCHX.EO.OlGO TO 4775 

JMX 4380 

00203* 


GO TO 6000 

JMX 4390 

00203* 

4775 

DPP 1 *DP2 

JMX 4400 

002040 


YB1*Y(N) 

JMX 4410 
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n r ? y 4 2 


r)p?3i)P2* • 9 

JMX 

4420 


740 

ICHX*2 

JMX 

4430 

n r ? 4 4 


GO To 6001 

JMX 

4440 

rr?T4^ 

f 000 

VS« ABS<YOUT-YRU 

JMX 

4450 



VNF x ABS ( YOUT -Y { N) } 

JMX 

4460 

OApn^l 


IF (VNE «GT • VS) GO TO 335 

JMX 

4470 

002*57 


Y«?*YH1 - -r. . . v ■ 

JMX 

44R0 

002**0 


OPPpaOPRl 

JMX 

4490 

OOP'i^l 


Yfl 1 =Y ( N J 

JMX 

4500 

002**3 


QPH1=DP2 

JMX 

4510 

00?^4 


GO T n 336 

JMX 

4520 

0O?O6 a 

?*»5 

YB?xY (N> 

JMX 

4530 

0*?*6* 


r>PP2=DP2 

JMX 

4540 

on?*7o 


0P2*0PBU ( YOUT-Y01) • (0PP3-DPBli / < YB3-Y81) 

JMX 

4550 

002101 

. 

GO TO 6001 

JMX 

4560 


C 

itfpation fnos 

JMX 

4570 


C 


JMX 

4580 

0(i? 1 o \ 

*73? 

IF<NSTEP.E0*1 >0Pl«0P3 

JMX 

4590 

0021 0* 


PCUM=PCUM+ ( (0P3*0Pl)*DX/?,0) 

JMX 

4600 

0021 1? 


JaNTT 

JMX 

4610 

Oo? 1 1 3 


RZ=R? U)*(U(3 t J#i )-U <?• J)) *<UC2 • J) -U (2 t J-l ) 1 

JMX 

4620 

C0?1?1 


R2sR7*tl (J1*<VIS< J) ♦ RHO ( 1 . J) *E ( J) i 

JMX 

4630 

0 0? 1 30 


UflRe ( Y C J ) »4 (NP*nPP)-Y (Ml **<NP*NPP) ) *DPI/ < (2 . ♦FLO* T (NP) ) #RHO 

JMX 

4640 



*(ltJ>4Y(N)44NL)-8Z/l«H0Mf J1*Y(N)**NL) 

JMX 

4650 

002161 


URRaUCLI* SOHTtURR) 

JMX 

4660 

0021 *4 


TOL 0 (2 ) * URR/UCtl 

JMX 

4670 

00216* 


O^LYaY (N) -Y {NN) 

JMX 

4680 

002170 


URfFxSQPT l(U(l ,nN)/0ELY-0ELY*DP1/4,/VIS<NN> )#RH0(1,NN)/ 

JMX 

4690 



•VIS (MN) ) 

JMX 

4700 


C - • 


JMX 

4710 


C . c f 

* * SECTION 04#* 

JMX 

4720 


C - 


* JMX 

4730 


c 

replace m line values with m.* line values 

JMX 

4740 


r ; 


JMX 

4750 


c * 

* • SECTION 9 # # # 

JMX 

4760 


c > 

# output section * * 

JMX 

4770 


c; .... 


JMX 

4780 

002203 

1 745 

XIN * xro*yt/(?.*ren) 

JMX 

4790 

002?0 7 


IF< NSTEP *EQ • 0 ) GO TO 4^8 

JMX 

4« 00 

002210 


IF(NSTEP .EG. 1) GO TO 908 

JMX 

4810 

00221? 


IF (XIN ,GE. PI1 GO TO 908 

JMX 

4820 

002214 


60 TO 7142 

JMX 

4830 

002215 

Q"8 

CONTINUE 

JMX 

4840 

002215 


IL a IL ♦! 

JMX 

4850 

002217 


Y T N» Y(N)/YT 

JMX 

4860 

002221 


UCL * U ( 2 » 1 1 *UCL I 

JMX 

4870 

002223 


TCL * T<?*1>*(TWREF-TCLI> ♦ TCLI* UCL # *2/ 1 2000 460. 

JMX 

4880 

002233 


AIJOM(I) aOf 

JMX 

4890 

002234 


oo nil i * 2 . n 

JMX 

49q0 

002235 


I M 1 « J-] 

JMX 

4910 

002237 


AUOM(I) « AUGMUM1) tCU<2»I)«UI2t!MlMMP5(I)«P§<lMin»(PSin- 

JMX 

4920 



1 PS(IMl)) 

JMX 

4930 

002251 

111! 

CONTINUE 

JMX 

4940 

002254 


Auom^t » AIJGM(M)/{P5NJ*2.*6,2836*#NL) 

JMX 

4950 

00226? 


PH?0 a PCtJM*VHFAD **1923 

JMX 

4960 


77 


an,, VISION — PSH LFVFL 29R-- NAS 


0*226* 


PFL ■ DFLTA*RNU/ (UCLI*OORFF) 

JMX 

4970 

t)n??7A 


W*!U<ft,ft3> It ? * IN f XRO » y 1 N * PH20* UCL * TCL t AUG^FT, URftl 

JMX 

49A0 

on?Mft 

ft3 

FORMAT {/t5X* UtfiKl3.3) 

JMX 

4990 

oc?iu 


PT ■ Pi ♦ #25 

JMX 

5000 

0ft?3?o 


DO 57 la 1 *MK 

JMX 

5010 

ro? - ^? 


IF { X I N .GE. PRoFEIXM Go TO 902 

JMX 

5020 

or?3?5 

57 

CONTINUE ‘ ~ ' ’ 

JMX 

5030 

002327 


GO TO 7142 

JMX 

5040 

002327 

4*8 

CONTINUE 1 

JM t 

5050 

002327 


WR I TE ( 6 . 752) 

JMX 

5060 

002333 


DO Qo9 J»1*N*1 — • ■ ’ 

JMX 

5070 

002335 


YPL(J) » YUJ/YOUT 

JMX 

5080 

00?337 


UJ( J) « U ( 2 ♦ J) *UCL I 4 ' 

JMX 

5090 

00234? 


LVH ■ UJ(J)**2/l2000.-460. * 

JMX 

5100 

002346 


TTP ( J) *T < 2 * J ) •(TWREF-TCLI) ♦ TCL 1 

JMX 

5110 

002353 


TTT ( J) *TTP < J) *LVH 

JMX 

5120 

00235ft 

909 

WRITE (6 i 61) j*YR L (j) »Uj ( j> *TTT(j) 

JMX 

5130 

002374 


WR T TE ( 6 # 757 1 

JMX 

5140 

002377 


60 TO 7142 

JMX 

5150 

002400 

902 

CONTINUE 

JMX 

5160 

002400 


PROFEU) a 50 • 

JMX 

5170 

Qft?40? 


WP I TE ( ft ♦ 351 ) 

JMX 

5180 

00240ft 

351 

FORMAT (//. 15X»*VEL0CITY AND TEMPERATURE DISTRIBUTION ARE PRINTED 

JMX 

5190 


1 

| OUT AT THE CORRESPONDING X IN. POSTION*t/) 

JMX 

5200 

00240ft 


wR I TE ( 6 * 752) 

JMX 

5210 

00241? 


00 709 J*1 » N 

JMX 

5220 

002414 


YPUJ) =» Y(J)/YOUT 

JMX 

5230 

00241ft 


UJ(J) *U ( 2 * J ) *UCL I 

JMX 

5240 

002421 


LVH = UJ(J) **2/12000.-460. 

JMX 

5250 

002425 


TTP ( J) *T (2 , J) * ( TWREF-TCL ! ) ♦TCL I 

JMX 

5260 

002432 


TTT U)*TTP( J) ♦LVH 

JMX 

5270 

002435 

7n9 

CONTINUE 

JMX 

5280 

00?437 


DO 1213 J* I tN» 2 

JMX 

5290 

002440 


WRITE (6,61) J * YRL ( J> »U J ( J) t TTT<J) 

JMX 

5300 

002453 


1= JM 

JMX 

5310 

002455 


iru .EO* N) GO TO 4311 

JMX 

5320 

002457 


GO TO 1213 

JMX 

5330 

002457 

433 1 

WRITE (6 >61) I»YRL(I) •UJ(X)«TTTU> 

JMX 

5340 

002473 

1213 

continue 

JMX 

5350 


r the 

FOLLOWING CARD is inserted IF a MASS *ENERGY CHECK REQUIRED 

JMX 

5360 


c 

CALL CHECK (YfE.U*RHOtYT.FOUCT.TTPiTTT#NN.UCLI *UJR) 

JMX 

5370 

00247ft 


CALL 0LChpK {N.U»R 5»Y,DX » DP2* UCL I ♦ RRH f TOLD • URR »RhO t AA ) 

JMX 

5380 

002511 


CALL CHECK (Y,E*U*RHOfYT,FDUCT»TTPfTTT»NN,UCLI,UJR) 

JMX 

5390 

002524 


WR I rE C 6 * 757) 

JMX 

5400 

002530 

7142 

DP 1 1 =DP 1 

JMX 

5410 

00253? 


DPI »DP2 

JMX 

5420 

002533 


DO 17 J*ltN 

JMX 

5430 

002534 


U < 1 »J)xU(2*J> 

JMX 

5440 

002540 

17 

TO t J) a T (2 * J) 

JMX 

5450 

002544 


TOLDU) * TOLD ( 2) 

JMX 

5460 


C 


JMX 

5470 


C * 

* • SECTION 10 * * * 

JMX 

5480 


C 

• TERMINATION TEST * * 

JMX 

5490 


C 


JMX 

5500 

002545 


IF (X ,GT. XX) GO TO IB 

JMX 

5510 


78 



RUN VrPSfrtN ?• 

i -- psr irveu 29*— 

NAS 




002551 


GO TO loo 



JMX 

5520 

C 




JMX 

5530 


r 

RETURN TO $T ART Of 

OUTER LOOP 


JMX 

5540 


c 




JMX 

5550 

002551 

l* 

CONTINUE 



JMX 

5560 

002551 


00 TO 8888 



JMX 

5570 

00255? 

883 

MR! Tf <6.??1 7) 

r AD u A T 1 >>>• fun 



JMX 

5580 


77 17 


\jr C ALCUL ■ T I UN w f f w 



55^0 


c 




JMX 

5600 


c 

format statements. 



JMX 

5610 


c 




JMX 

56? 0 

002555 

61 

FORMAT <25x, I6.6X.F 10.6*6x t Fl4*5.2KfF14.5) 


JMX 

5630 

002556 

767 

FORMAT {// *6X * * I 

X IN. X/H 0 

B/BO 

P JMX 

5640 


•* 

1 h20 UCEnT(F/S> 

TOCEnt (OEij.R) augment 

USTAR(F/S)*./) 

JMX 

5650 

002556 

752 

FORMAT ( 2S* *• 

J V(J> 

Ut J) 

JMX 

5660 



1 TIJ) •> 



JMX 

5670 

002556 


STOP 



JMX 

5680 

002560 


£Nf> 



JMX 

5690 



«UN 


yrRSIO*' ?»3 --PSR LEVEL 298-- 


naS 

program length including; I/O RdFFERS 

01341 * 

FUNCTION ASSIGNMENTS 


F 

- 

000006 




STATEMENT 

ASSIGNMENTS 



7 

- 

00264s 

9 


000500 

12 

. 

0012*3 

1 6 


001742 

22 


O00241 

23 


000266 

39 


P02635 

45 


001044 

63 

- 

003013 

73 

- 

002643 

117 

. 

000444 

118 


000421 

219 

• 

001477 

235 


002064 

351 

- 

003020 

458 


002327 

740 

- 

002043 

752 


003067 

778 

• 

001047 

796 


001011 

83P 

. 

000035 

883 


0O2552 

1 02 1 

- 

000054 

1051 


002744 

1360 

- 

001007 

1373 


000760 

2001 

. 

000211 

4172 


001076 

4775 

- 

002036 

6000 


002045 

714? 

• 

00?53o 

7232 


002752 

7817 

- 

002631 

8019 


002577 


BLOCK NAMES AND LENGTHS 
- 000004 


Variable alignments 
A - poGHTS AA 


005621 

AMR 

- 

007221 

ASTAR 


007175 

AZ 

• 

007 , 0 S 

AZ1 


007377 

BB 

• 

007312 

BB 8 


006465 

BH 

• 

nr*7234 

BLEND 


007257 

B7 

- 

00*341 

C 


006251 

CON A 

• 

007265 

CONfl 


007266 

CONI 

- 

007?7i 

D 


oo57?7 

OELTA 

- 

007242 

DELU 


007302 

Djr 


007272 

DPB 1 


007333 

DPU 

- 

007156 

OP? 


007155 

El 

. 

004541 

FDUCT 


007245 

HALF 

• 

007204 

I 


007162 

IFLOW 

• 

007237 

IL 


007230 

IP1 

• 

007300 

IT 


007304 

J 

. 

007252 

JFLOW 


007231 

JP1 

• 

007253 

KJ 


007243 

mk 

• 

007161 

N 


007222 

NN 

• 

007224 

NP 


OOOOOOCOl 

ntest 

- 

007161 

NTK 


007315 

PA MR 

. 

007146 

PAi 


007210 

PE 

- 

00*217 

PH20 


007235 


NAS 


{ 


10 


001216 

11 


OOU65 

lfl 


002551 

19 

" m 

000262 

29 


000256 

37 


002633 

48 


001221 - 

6 ] 


003044 

88 


0007U 

100 


000670 

119 


001447 

162 


000726 

236 


002070 

279 


002603 

631 


001221 

737 


000045 

757 


003051 

762 


001013 

629 


002616 

832 


001163 

902 


002400 

905 


002215 

1201 


000216 

1213 


002473 

1748 


002203 

1934 


001527 

4173 


001116 

4311 


002457 

6001 


001330 

673? 


002101 

7347 


001305 

7717 


003035 

8140 


001132 

6888 


000016 


AMaSSO 


007171 

AMASS! 

- 

007170 

AUGM 


006575 

AUGmET 

- 

007347 

A p 


OO 7 I 74 

B 

- 

O06W3 

R DELTA 


007275 

BE 

- 

007274 

ROUND 


007330 

BY 

- 

00 * 22 * 

CC 


00*153 

CHEki 

- 

00*205 

COND 


00*270 

CONE 

- 

00*26* 

obar 


007327 

OEL 

- 

007350 

DFLY 


007342 

DIG 

- 

00*326 

DPP2 


007340 

DPI 


007154 

OX 


007311 

E 

- 

005063 

GAmA 


007146 

H 

w 

006357 

ICHX 


007314 

ICORE 

- 

007241 

I M 1 


007301 

INM 

- 

007276 

ITER 


007226 

IUftfS 

• 

007233 

JJK 


007273 

JM 1 

- 

007254 

LVH 


007351 

LZ 

- 

003244 

NJ 


007251 

NL 

- 

000003C01 

NPP 


000002C01 

NSTfP 

- 

007236 

NTP 


007256 

NTT 

- 

007223 

PA 2 


007214 

PCUM 

• 

007264 

PT 


007225 

POl 

- 

007164 


80 



run VFR510N 2.1 — PSR LEVEL 2RA— 


pp 

- 007)4? 

PREF 

007145 

PS 

• 0043?* 

PSN 

0072*6 

rd 

- 007172 

REN 

0072** 

RM? 

007206 

RNU 

0071*3 

R2 

004413 

SUP 

000001C01 

S3 

004217 

T 

(•03567 

TESTl 

007202 

TEST? 

007203 

tom 

006*73 

TOO 

007167 

r sec 

00721* 

TTP 

005513 

U 

003245 

UCL 

0073*5 

UJR 

007220 

UN 

007307 

UPR 

007240 

USEC 

007216 

VNE 

007336 

vs 

007335 

X in 

■' 007343 

XRO 

007261 

X22 

007200 

Y 

00*755 

YD 

007310 

YIN 

0073** 

YPPT 

007332 

YRL 

007035 

YTEST 

007331 

Y 1 

007176 

Y4 

Yfl 

007320 

007324 

Y5 

Y9 

007321 

067325 


START OF CONSTANTS 
002562 


start of temporaries 

003100 

START OF IwOIRECTS 
003166 


NAS 


proPe - 

006703 

PRT 

m 

007150 

PSNJ . 

007247 

PI 

m 

007207 

RHO 

00^171 

Rm 

m 

007211 

RORCf - 

007144 

Rl 

m 

004647 

si 

004003 

S2 

m 

004U1 

tct 

007346 

TCLI 

m 

007212 

t JR 

00^263 

TOL 

m 

007201 

T01 « 

007165 

TREF 

m 

007151 

TTT 

003461 

TWRfF 

m 

007152 

UCL! - 

007213 

UJ 

m 

006727 

UNN 

007306 

UftEF 

w 

007232 

VHC40 7 - 

007262 

VIS 


005405 

X 

007157 

X8LENQ 


007260 

XX 

007160 

xu 

m 

007177 

Y01 

007334 

YB2 

m 

007337 

Y J 

007173 

YOUT 

m 

007256 

VS 

007303 

YT 

m 

007255 

Y2 

007316 

Y3 

m 

007317 

Y6 

Zk * . 

007322 

007313 

Y7 

m 

007323 


UNUSEO COMPILER SPACE 
062200 



pus version 


2 .3 — PS« LEVEL Z9P 


000" I * 
ooo* i o 
000 n l 0 
000*1 1 
OOOOJ 7 
oof>n?3 
00004* 
000050 
000050 


500 


SUBROUTINE Yl)iS(Y f pS*HHntU<Ni 

Vis 

10 

DIMENSION Y(7o) ,PS(70) . PhO < 2, To) *U <*, ?0) 

YIS 

30 

COMMON NP»SOK*NPP t NL 

YIS 

40 

Yd > *0.0 

YfS 

50 

no 5o o i«?»n 

YIS 

60 

Zs(?.*ELOAT(NPn/(RHO(?f 1-1) *U<2'!*W ♦«M0 <?♦ 11 # U<2* Ill 

YIS 

70 

7y»y( i-ii*#<Np*Npp) ♦ z«<p$d>*psti>-psu-n # p$u-in 

YIS 

80 

V(T) a Z Y**SOP 

YIS 

90 

RETURN 

YIS 

100 

END 

YIS 

HO 


RUN VERSION Z.3 --P5H LEVEL 29A*« 


YDl $ 


ms 


SUBPROGRAM LCNATH 
0001 0* 

FUNCTION ASSIGNMENTS 

STATEMENT ASSIGNMENTS 

BLOCK NAMES ANn LENGTHS 
- 000004 

VARIABLE ASSIGNMENTS 
I - OftOl 03 NP 

z - 000104 ZY 

START OF CONSTANTS 
00005? 

START OF TFmPOBaHIES 
000054 

START OF TN0IRECT5 
000073 

UNUSED COMPILER SPACE 
073300 


oOOOOOCOI NPP « 000002C01 SOP « OOOOOICO 1 

000105 


82 



HU* version 2.3 — PSH LEVFl 29*-* 


oooon 
0000*1 
00001 ? 
0000*4 
oooon 
000017 
000020 
000021 
00002S 
000031 
00C 043 
000044 
00004* 
000050 
000051 
00006? 
000062 


C 

C 

C 

C 


1 


? 


SUBROUTINE CAU <A »B'<; tft 


45 P** MAY *9 |9Tl 


* ♦ 1 THIS EVALUATES *ESULT$ U5J*0 THOMAS ALGORITHM ♦ 

0ImEn$Jqn *(70 ) ?B (70 ) i C ( 70 ) tO ( 70 )*H(70 >»*(70 )*O(70 )»G(70 
N 2-U-2 
Nl*J-l 


win.AU) 

G ( 1 ) «D ( 1 ) /W ( * I 

no 1 k* 2 *n? 

K 1 »K •• * 

0(Kl ) *B ( K 1 )/W(Kl ) 

y(KJ,A(K).C(K)#0(K*) 

G(K)p(0(K)«Q(K)«Q(K*))/W(K) 

H(m2)«G(N2) 

N3rJ-3 


00 2 K.* t N3 
KK»N?-K 

M(KK)»G(KK)-Q(KK) *HCKK*t I 
RETURN 

end 


CAL 

CAL 

CAL 

CAL 

cal 
) CAL 
CAL 


CAL 

CAL 

CAL 


ca l 

CAL 

CAL 

CAL 

cal 

CAL 

cal 

cal 

CAL 

cal 

CAL 

CAL 


10 

20 

30 

40 

50 

60 

70 

00 

90 

100 

no 

120 

130 

140 

150 

160 

170 

100 

190 

200 

210 

220 


n\Jh' VrPSIOK) J> # 3 — P$» level, 39 A-* 


CALC 


CALC 

SUBPROGRAM LENGTH 
000437 

r uNCTtoN assignments 

STATEMENT ASSIGNMENTS 

8lOC* NAM£s ANH LENGTHS 

VARIABLE ASSIGNMENTS 
G * OftO??3 K 

N 1 - 000432 m2 

V - 000107 

START OP CONSTANTS 

000064 

start of temporaries 
00006 $ 

START OF 1NMRECTS 
000071 


000433 KK * 900436 K* 

00043* n3 w Q004J5 9 


000434 

000215 


UNUSen COMPTIER SPACE 
07320ft 



RUN VFosION ?,** --P5W LEVEL 298 — 



- 

SUOPOUTINE PROF<U # T f NJ,N*UJRf TJR,KJ,i 

UCLI) 

PRO 

10 


C 

VERSION FOP SINGLE STRfaM FLOW WITH 

a wall boundary layer 

IF K J* IPRO 

20 


c 

VERSION EO» TWO STREAM FLOW (TOP-HaT 

PROFILE ) OR SINGLE 

PRO 

30 


c 

STREAM FLOW wITmOUT A BOUNDARY Layer 

IF Kj«2 

PRO 

40 


c 



PRO 

50 

000013 


n I MENS I ON UI2t7f>)tT<2?70> 


PRO 

60 

ooonn 


IF ( K J * EO * 1 ) GO TO 12 


PRO 

70 

ononis 


nO lo I ■ 1 * n 


PRO 

80 

000016 


TFU.GT.NJ ) GO TO 5 


PRO 

90 

000021 


l)(M)*1.0 


PRO 

100 

000023 


TO *T)»0.0 


PRO 

110 

000024 


GO TO 10 


PRO 

120 

000025 

5 

U«1 tT)»UJR 


PRO 

130 

000030 


7(1 i T)«TJR 


PRO 

140 

000031 

10 

CONTINUE 


PRO 

150 

000034 


GO TO 9 


PRO 

160 

000034 

12 

REAP (5 f 1 3) (U(l f 1) f I«l f N> 


PRO 

170 

00005* 

13 

FORMAT (6E13.6) 


PRO 

180 

000054 


DO 15 I *1 * N 


5 PRO 

190 

00006) 

15 

T(i»n«o*o 


PRO 

200 

000065 

8 

nO 2o 1*1 *N 


PRO 

210 

000067 


U (2 * I ) »U (1*1) 


PRO 

220 

000073 

20 

T<?*T>*T<l»I> 


PRO 

230 

000077 


U(]|N)*0.0 


PRO 

240 

000100 


IJ ( 2 1 N) »0 * 0 


PRO 

250 

000102 


RETURN 


PRO 

260 

000102 


END 


PRO 

270 

RUN VfpSTAN 2.3 

— PSR LEVEL 298 — 

PROF 




PROF 

SURP p OGP^M LENGTH 
COM ?4 

FUNCTION ASSIGNMENTS 

STATEMENT ASSIGNMENTS 
5 - *00025 6 

13 - enoin? 15 

block names an o lengths 

Vac?TA9L f assignments 
I - 0001*3 KJ 

Start of constants 
000104 

START OF TEMPORARIES 

coot n 

START of IMOIRECTS 
000 H f 

UNUSEO CCMPTtfO SPACE 
073100 


000065 10 - QQ0031 12 - 00003* 

000061 


oooooo UCLI * 000001 


84 



3 U *; VERSION ?,1 — PSR LEVFL 2 98 — 


SUp&oUTINE PSI( N t PSN.NJ.PS. PSNJ) 


c 


00001 " 


niMFNSION PS ( 70 1 


C 

JN *1 IS THE NUMBER OF NODE POINTS IN JET 

OOOOl o 


JN * 11 

00001 o 


NJ* JN* 1 

OOOOl ? 


psm *n.o 

000013 


DO 30 I * ?*NJ 

000015 


PS (I ) * PSnJ/2.*<1 ,-005 <3.1416/FloAT(Jn) # Pi 0*T(I-1 1 ) ) 

00003? 

30 

CONTINUE 

000034 


FS ■ 0.1*<PSN-P5NJ) 

000036 


J «N J* 1 

000040 


K* N J ♦ JN 

00004? 


no 40 I * JfK 

000044 


PSU) »PS( NJ! ♦ FS*(1. -COSO. 1416/(2. *JN)#FLOAt(I-NJ))> 

00006? 

40 

CONTINUE 

000065 


OFLPS *.85*(PSN-PSNJ)/FL0AT (N-K-181 

000071 


J *K ♦ 1 

000073 

'■ '■ 

K * N-18 

000075 

' ■ - 

no 50 I « J »K 

00007? 


PSU) *PS(I-1) *0ELPS 

00010? 

50 

CONTINUE 

000105 


FSA*.05*(PSN-PSNJ) 

000107 


J * K + l 

000111 


RFAO (5.200) (PSU) tI»J*N) 

000130 

?00 

FORMAT (6E1 3.6) 

000130 


DO 60 I » J t N 

000135 


psu) • ps(k) ♦ FSA*Ps<n 

00014! 

60 

continue 

000143 


DO 101 1*1. N 

000145 


PS<t ) * SORT Cps C i ) > 

000155 

1 Ol 

CONTINUE 

000160 


RETURN 

00016" 


END 

RUN VFOSTCN 2 » 

3 — P$R LEVEL 298— PSI 


PS! 


PSI 

PSI 

PSI 

PS! 

PSI 

PSI 

PSI 

PSI 

PSI 

PSI 

PSI 

PSI 

PSI 

PSI 

PSI 

PSI 

PSI 

PSI 

PSI 

PSI 

PS! 

PS! 

PSI 

PSI 

PSI 

PSI 

PSI 

PST 

PSI 

PSi 

PSI 

PSI 

PSI 


SUBPROGRAM LENGTH 
000?26 


FLECTION ASSIGNMENTS 


S T fl r E k *FNT ASSIGNMENTS 
200 - 000174 

BLOCK N4V£s ano LENGTHS 

V4PTio' ta e ASSIGNMENTS 

OFL 3 S - Ono??4 FS * C00221 F54 » 0002?5 I - 000220 

J - PO02?2 JN - 000217 K • 000223 

Stact of CONSTANTS 
00016 ? 

START OF TpvPOOaPIES 
00 0 1 76 

Start of InhIRFcTS 
0C0"1 s 


10 

20 

30 

40 

50 

60 

70 

80 

90 

100 

no 

120 

130 

140 

150 

160 

170 

180 

190 

200 

210 

220 

230 

240 

250 

260 

270 

280 

290 

300 

310 

320 

330 


unused compiles space 

07270" 



RUN VFRSIDn ?.3 — PSR LEVEL ?9fl 




SUBROUTINE Radius (X, tout tK7*RFN) 


GEM 

10 


C 

AS PER M*V 19 1971 


GEM 

?o 

000007 


DIMFNSION RR ( 25 1 ♦ XX (25) 


GEM 

30 

000007 


Go To 1 1 * 2 ) *K2 


GEM 

40 

00001* 

1 

READ ( 5 # 6 ) NS 


GEM 

50 

000022 

6 

FOpMAT (12) 


GEM 

60 

00002? 


READ ( 5 »3 ) ( RM < 1 1 f I ■ 1 * N$ ) 


GEM 

70 

000*35 


READ (5*3) < X X ( I ) ♦ 1«1 «NS) 


GEM 

80 

oonnso 

3 

FORMAT (8F1 0,0) 


GEM 

90 

000050 


WRITE (6,20) 


GEM 

100 

000054 

'?0 

FORMAT (//*?5X,*0UCT GEOMETRY* f // f 5X*#0UcT DUMFTFR 

OR HEIGHT 

GEM 

no 


1 DISTANCE From N0Z2LE EXIT (IN)*.//) 


GEM 

120 

000054 


DO 25 I* 1»NS 


GEM 

130 


C . 

the following caro is insert if half width qf the 

DU$T is gSEO 

GEM 

uo 

000060 


RR ( T ) *RR { I ) *? , 


GEM 

160 

OOOOf.? 


WO T TE ( 6 • 7 ) WHdlfXXm 


GEM 

1 70 

00007?* 

7 

FORMAT UOX«Flo 9 5«2SX 9 Fl 0.5) 


GEM 

180 

0000?? 

?5 

CONTINUE 


GEM 

190 

000077 


DO 4 I * 1 1 N5 


GEM 

2QQ 

000101 


RR(I)*RR(I)*REM 


GEM 

210 

000103 

4 

XX<i>»XX<i>*R£n*2. 


GEM 

220 

000107 


T B 1 


GEM 

230 

000110 

? 

T F(X.GT # XX(NS) )G0 to 5 


GEm 

240 

00011 A 


TFlX,GT.XX<l*i)>I»I*l 


GEM 

250 

000120 


YOUT^RH ( ! ) ♦ ( X-X* ( I ) 1 • (BB ( I ♦ 1 ) -RR ( l ) > / (XX ( I ♦ l ) -XX 1 1 ) ) 

GEM 

260 

000131 

5 

RETURN 


GEM 

270 

00013? 


END 


GEM 

280 


RUN VERSION ?.3 --P$R LEVEL 298-- 


RAOIUS 


RADIUS 

SUHPPOGPAM LENGTH 
OO0?66 


FUNCTION ASSIGNMENTS 


STATEMENT ASSIGNMENTS 
1 - 000014 2 

** 0001 10 

3 

- 000140 

5 

~ 000131 

6 <*> 000136 7 

- 000160 

20 

- 000143 



BLOCK NAMES AND LFNGTHS 






variable assignments 

I - 000265 NS 

- 000264 

RR 

* 000202 

XX 

- 000233 


START OF CONSTANTS 
OOOl^A 

START OF TEMPORARIES 
000164 

START OF TNOIRFCTS 
000174 

UNUSE" COMPILER SPACE 
073000 
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RUN VFPSTOm ?,3 — PSR LEVFL 29ft 



SUBROUTINE Check <Y,E •U*RHO*YT,FOUCTtTTPtTTT»NN,UCLl iU JR 1 

CHE 

10 

000O1 A 

common np,sup,npp,nl . . ■ - l 

CHE 

30 

000016 

DIMENSION Y (70) , E (70) *U(2»T0) f PHO(2 # t01 * f tP ( 70 1 i XTT ( 70) 

CHE 

40 

00 001 * 

ENFPr,*o. 

CHE 

50 

000017 

TFLOW* 0 , 

CHE 

60 

000020 

V I SC I N«0 • 

CHE 

70 

000021 

nEL? 3 0. 

CHE 

80 

00002? 

DO 34?1 J»1 iNN 

CHE 

90 

000027 

JPlaJ^I 

CHE 

100 

000025 

TFLOW «TFLOW ♦ f Y ( JP l ) ** (NP*NPP) - Y ( J) ** <NP*NPP) ) 

CHE 

no 


1 *(RhO(2* jPI ) ♦Rh'' (2* J) ) 

Ch e 

120 


?*(Mf? f JPlj *U<2, J n/ 4« 

CHE 

130 

000056 

VI5 C lN=VIScINMY(jPn*Y(J) )•(£ C JPl ) ♦£ ( J> > * (U (2* JPl ) -U (2« j) ) 

CHE 

140 


1 / ( T T P ( J ♦ 1 ) ♦TTP(J) )/2. 

che 

150 

0001 Op 

OEL?=DEL2* (U(2, JPl ) ♦U<2t J) ) # (V ( JPl) ** <NP*NPP> -Y ( Jl *4 (NP^NPP) ) 

CHE 

160 


1 * ( (U (2* JPl ) ♦U (2. J) ) * ( .5/u JR) *1 . ) *.5/UJ« 

CHE 

170 

000134 

FNFPG=ENERG* CTTT(J) ♦ TTT( JPl ) >* (Y ( JPl ) **(NP*NPP> -y( J) ** (NP*NPP) ) 

che 

lftO 


l*(ll(?t JPl ) ♦U(2,J> ) * (RhO (?, JPl ) *RH0 

CHE 

190 

00017? 

3471 CONTINUE 

CHE 

200 

00017*; 

IF (NP .EG* 0) r, 0 TO 219 

CHE 

210 

000176 

FNFRGsENERG/ (FOUCT*YT**?> 

CHE 

220 

000200 

tFlOwsTFLOw/ (F0UCT*YT**?1 

CHE 

230 

000201 

GO TO 319 

CHE 

240 

00020? 

? 1 9 CONTINUE 

CHE 

250 

00020? 

PSN a FOUCT ' y 

CHE 

260 

000203 

FNFPGaENERG/PSN y i? 

CHE 

270 

000205 

TFLOWaTFLOW/PSN 

CHE 

280 

000206 

319 CONTINUE 

che 

290 

00020* 

V I SC I N*V I SC IN*. 00 0158**? 

CHE 

300 

000210 

nEL?«DEL2 **.5*.000158/UCLI 

CHE 

310 

000210 

WHlTF(6>347?)ENE^G ( iFL0w*vJs^^Hf0EL2 

CHE 

320 

00023? 

147? FURMAT (/,6X**T0-F « * »Fi 4 .6 *6# i *TOTAL FLOW ■ *#El4*6» 

che 

330 


16 x ,* v ISCIN = *,El4 # 6,2x** n ^L?.*,El4 i 6) 

CHE 

340 

000?3? 

RETURN 

che 

350 

000233 

FND 

CHE 

36o 


RUN VPOSTOM ?.3 — P$R LEVEL 29ft — 


check 


CHECK 

SijppopGOAM LENGTH 
000409 


FUNCTTON ASS I GNmENTS 
5 T ATE V FNT ASSIGNMENTS 


219 

0002"2 

319 - - 

000206 

3472 

- 000245 



8L0CK NamfS AND LENGTHS 



- - — 

■ 

- 

- 

0 0 0 0 04 







VARJ A9LE 

ASSIGNMENTS 






CEL2 - 

000401 

ENERG - 

000376 

J 

- 000402 

JPl 

- 000403 

NN 

000002 

NP 

oooooocoi 

NPP 

- 000002CO 1 

PSN 

- 000404 

TFtnw . 

000377 

TT P 

000000 

TTT 

- 000001 

UCLT 

- 000003 

UJP 

000004 

VISCIN - 

000400 





ST A PT OF 
C C 0235 

CONSTANTS 







START OF 

TEMPORARIES 






00026" 








START OF 
0003^4 

INOIPFcTS 








UN'JSFn COMPILE® SPACE 
072500 



RUN VPDSTOn 2*3 — PS« LEVEL 29B-- 


000017 

000017 

000017 

000020 


00003* 

000037 

0000*1 

00004? 

1 

OOOO** 

000061 


00007S 
0001 00 
00011* 
000121 
00012* 

2 

0001*? 

*0 

OOOi *2 
0001*3 

i 00 


StmROUTlNF WLCHFK (N,UfPS,y # OX *0P?*UCLI fH 
OT^FNSIDN UI2t 7 0) tPS< 70) • DUX < 7O ) t RHO < 2 , 7 0 ) 

* tTHL0{2> 

COMMON NP • SGP « NPP »NL 
no 1 J ■ 1 *N 

nux(j) m (U( 2 tJ)-uu»J))/ox 

CONTINUE 

A & ( 1 ) a 0 • 

no ? j« 2 * n 

JM1 * J-l 

A A ( J) «AA < j-l ) ♦ ( nux ( J) *PS ( J) ♦OUX ( JM] } «P5 ( JMj ) 
RRPU) =»0P? *y U> •* (NP + NPPJ / (2,*FL0AT (NP + N 

CONT I NUE 

r = -PhO <2,1 ) * ( TOLD (2 ) ♦ TOLOU) ) **2*Y( 

AORsRBB (N) / C ti - A A (N) ) 

0C * B-AA(N) 

WRITE ( 6* *n ) BC * RBB (M • ADR 

FOPm a T </* l 0 X,#fLUId FORCES * **£13.6** 

1 Eli. 6** RATIO OF THESE TWO • **E13*6* 

RETURN 
END 


left* Told * urr * rho * a a i 

BLC 

10 

*AA (70) *y < 70 J ffiPH C 7 0 ) bLe 

30 


BLC 

*0 


BLC 

50 


BLC 

60 


BLC 

70 


BLC 

BO 


BLC 

90 


BLC 

100 


BLC 

ll(t 

)*(Psui-PS(jMm 

BLC 

120 

ippn 

BLC 

130 


BLC 

1*0 

N)»*NL«.25 

BLC 

150 


BLC 

160 


BLC 

170 


BLC 

1 BO 

press, forces • •* 

BLC 

190 

/) 

BLC 

200 


BLC 

210 


BLC 

22 0 


HU f ‘ 


VFQSTAN 2,1 — PSR LEVEL 29R-- 


BLCHEK 


RLCHEK 

SUBPROGRAM LENGTH 
0003*2 

FUNCTION ASSIGNMENTS 

STATEMENT assignments 

*0 - nooiM 100 

BLOCK NAMES ANO LFNGTHS 
- .000004 

VAPTAO(_F ASSIGNMENTS 
AA - 000005 AOB 

BC - 0001*1 DUX 

NL - 000001C01 NP 

TOLD - 000002 UCLI 

START OF CONSTANTS 
0001*S 

START OF TEMPORARIES 
000165 

START OF TNDIRFCTS 
000?]o 

unused compiler space 
072700 


0001*2 


0003*0 R 

000227 J 

OOOOOOCOl NPP 
000000 URR 


000337 BBB - 000001 
000335 UM1 * 000336 
000002C01 RHO * 00000* 
000003 
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RUN VrRSTOw 2,3 --PSR LEVFL 29fl 


OOOO ?o 

0000?« 

0*0n?1 

000023 

00002* 

000044 

00004* 

000050 

00005? 

00005* 

000057 

0000*? 

0000*5 

000070 
00007? 
000073 
000103 
00010* 
000110 
000113 
00011 4 


000133 
0001 3* 
000 1 * n 
000144 

ooou* 

00016* 

000171 

000173 

000177 

000 ? 0 ? 
000203 
000?11 
000?1 4 
000?1* 
000??1 

000?3* 

000?40 

000243 

00024* 

O00 ?47 


SUPPOUTINE LOOK | J JK ,NN, U, Y,OEtT A, BH, BY, N,CC.lFLOw.XBLEND,Yj, BE) 

C FOR F JECTOP FLOW AS PER FEB. 19 

C 

DIMENSION U (g *70) t Y 1 70 ) 

C S ea rCH F 0 p InFLFCtIo n Pn^T 

no 770 J* J JK *NN 
M* J ♦ t 
KaJ*? 

nZ«(UU »J)-U(1*M) )/<Y(J)-Y(M) )-(U(l»M)-U(lfK) >/<Y(M)-Y(K) 1 
IF(OZ,fiT.O.)GO TO 771 
770 CONTINUE 
771 IFI07.LT. o*>60 TO 776 
c no inflection point 

jj«n-m -1 
llPOT aU ( 1 * M) 

. TFtY(M) ,LE.BY)0O TO 776 
OE!.TA*Y (N) -DELTA 
TF(Y(M),6E, BELT A ) 00 TO 776 
C SEARCH FOR SHEAR LAYEP OUTER EDGE 

mY»m-1 

no 1319 1*1 * MY 

UTY*<U<1 tM-I)-UPOT>/(UM ,11-UPOT) 

TFfUTY.r,T..003)GO TO 13?0 

1319 CONTINUE 

1320 MP*M- I ♦ 1 

JJK*MP 

RY* < Y (MP-1 ) -T ( Mp) ) * ( * 00 3* ( U ( 1 1 1 1 -UPOT 1 * 

UMl *MP-1> *UPOT)/ (UO fHP-1 )-U(l ,HP) ) ♦ Y (MP-1 ) 
c SEaPcH FOR Shear LAYER INNER eO&E 

00 77? I *3 , M 

UTT*(UU *l)-UUf t) )/(UM tl)-UP0T) 

TFfUTT.OT*. 003)00 TO Th 
77? CONTINUE 

773 re*ym-i) ♦(yid-vih) >*< un 

1 s oo3»(uut n*uPOT>)/(UMt D.unfi.n) 

IFU.E0.3)RE»0. 

RFPAC*KE/Yj 

IF (PFPAC.LE. #050) CC* • 1 0* 

F1H*RY-BE 

C SEARCH FOR EDGE OF WALL BOUNDARY LAYER 

00 774 1*1, JJ 
UKK*<UPOT-UU#M*|) J /UPOT 
TF( UK K.fiT., 010) GO TO 77? 

774 CONTINUE 
775 Kl*M*I-i 

RE) TAsY (KI) -<Y(KI«1)-Y(KT) )• IU(1 'KI ) -»990« 

1UP0T)/(U(1 ,KI*1) -U(l fKT) 1 
IF(RV.GE.PELTA)GO TO 77* 
nELTA*Y (N) #*8ELTA 
MZaK f -HP 

IF (H7.LE.1 )GO TO 776 
r, 0 TO 777 

C TFLOwl THERE IS NO INFLECTION P0IN T 


LEF 

10 

LEF 

20 

LEF 

30 

LEF 

40 

LEF 

60 

LEF 

70 

LEF 

BO 

LEF 

90 

LEF 

100 

LEF 

no 

LEF 

120 

LEF 

130 

LEF 

140 

LEF 

150 

LEF 

160 

LEF 

170 

LEF 

180 

LEF 

190 

LEF 

200 

LEF 

210 

LEF 

220 

LEF 

230 

LEF 

240 

LEF 

250 

LEF 

260 

LEF 

270 

LEF 

2A0 

LEF 

29 0 

LEF 

300 

LEF 

310 

LEF 

320 

LEF 

330 

LEF 

340 

LEF 

35o 

LEF 

3*0 

LEF 

370 

LEF 

380 

LEF 

390 

LEF 

400 

LEF 

4)0 

LEF 

420 

LEF 

430 

LEF 

440 

LEF 

450 

LEF 

460 

LEF 

470 

LEF 

4R0 

LEF 

490 

LEF 

500 

LEF 

510 

LEF 

520 

LEF 

530 

LEF 

540 

LEF 

550 


RUN VFPSTON 2.3 —PSR LEVEL 298— ~ COOK 


00025" 
C 0025? 
0002*3 
00C254 


776 TFL0W*1 
xBLFND*0* 
777 PETURN 
END 


LEF 560 
LEF 570 
LEF 580 
LEF 590 



RUN VFBSTON 2.1 L€VFL 2o«--» LOOK 


UOOK 

SUBPROGRAM LENGTH 

00012s 

FUNCTION ASSIGNMENTS 
STATFmFNT AS$ T ftNMFNTS 


771 

0000^0 

773 

w 

000146 

775 

• 

000216 

776 


000250 

777 

ono2Si 

1320 

- 

ooono 







BLOCK NamFS ANO LENGTHS 









variable 

assignments 









be 

nooooAj 

GELTA 

m 

000313 

0FPAC 

• 

000321 

BY 

- 

000000 

cc 

000002 

□ 2 

m 

000310 

I 

m 

000315 

IFLO* 

„ 

000003 

J 

oooins 

J J 

s> 

O003U 

K 

m 

00030? 

K I 

• 

000323 

M 

000304 

HP 

* 

000317 

MY 

m 

000314 

MZ 

» 

000324 

N - 

0 oooo i 

UKK 


000322 

UPOT 

w 

000312 

UTT 

m 

000320 

UTY 

000116 

xbleno 


000004 

YJ 

m 

ooooos 





START OF CONSTANTS 
000256 

START OF TFMPO»aRIES 
000264 

START OF TN0IRFCT5 
000277 

UNUSFO C0MOTLFR SPACE 
072400 
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Table 1 


Mixing Section Dimensions for 1.875" Throat Size 


X 

Inches 

Inches 

X 

Inches 

±y 

Inches 

X 

Inches 

+7 

Inches 

-3.000 


-0.937 

1.960 

1.062 

1.280 

-2.937 

4.496 

-0.875 

1.922 

1.125 

1.268 

-2.875 

4.288 

-0.812 

1.886 

1.187 

1.256 

-2.812 

4.128 

-0.750 

1.850 

1.250 

1.244 

-2. 750 

3.982 

-0.687 

1.814 

1.312 

1.234 

-2.687 

3.860 

-0.625 

1.782 

1.375 

1.222 

-2.625 

3. 740 

-0.562 

1.746 

1.437 

1.210 

-2. 562 

3.628 

-0.500 

1.715 

1.500 

1.200 

-2.500 

3. 510 

-0.437 

1.688 

1.562 

1.192 

-2.437 

3.412 

-0.375 

1.658 

1.625 

1.182 

-2.375 

3.310 

-0.312 

1.632 

1,687 

1.174 

-2.312 

3.226 

-0.250 

1.608 

1.750 

1.162 

-2.250 

3.140 

-0.187 

1.586 

1,812 

1.154 

-2.187 

3.064 

-0.125 

1.562 

1.875 

1.146 

-2.125 

2.998 

-0.062 

1.540 

1.937 

1.138 

-2.062 

2.918 

0.000* 

1.519 

2.000 

1.130 

-2.000 

2.840 

0.062 

1.502 

2.125 

1.116 

-1.937 

2. 778 

0.125 

1.486 

2.250 

1.100 

-1.875 

2.708 

0.187 

1.470 

2.375 

1.086 

-1.812 

2.646 

0.250 

1.452 

2.500 

1.077 

-1.750 

2.586 

0.312 

1.438 

2.625 

1.066 

-1.687 

2.520 

0,375 

1.422 

2.750 

1.058 

-1.625 

2.464 

0.437 

1.406 

2.875 

1.048 

-1. 562 

2.404 

0.500 

1.394 

3.000 

1.042 

-1.500 

2.350 

0. 562 

1.380 

3.125 

1.038 

-1.437 

2.296 

0.625 

1.366 

3.250 

1.034 

-1.375 

2.246 

0.687 

1.354 

3.375 

1.030 

-1.312 

2.196 

0.750 

1.340 

3.500 

1.024 

-1.250 

2.152 

0.812 

1.328 

3.625 

1.018 

-1.187 

2.110 

0.875 

1.318 

3.750 

1.014 

-1. 125 

2.070 

0.937 

1.304 

3.875 

1.010 

-1.062 

2.030 

1.000 

1.288 

4.000 

1.007 

-1.000 

1.990 



8.000 

0.938 





11.000 

0.938 





23.500 

1.593 


T 

Nozzle discharge plane at x = 0. 000 
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Table 2 


Variation of Individual Integrated 
Traverse Mass Flows For Each Test Run 


Run 

Variation of Integrated Traverse 

Number of Traverse 


Mass Flow Rate Around An 

Locations 


Average Value 


1 

+2.2%, -3.0% 

4 

2 

+3.5%, -2.3% 

5 

3 

+3.4%, -3.9% 

9 

5 

+0%, -0% 

3 

6 

+2.3%, -2.8% 

4 

7 

+41%, -2.8% 

5 

9 

+5.6%, -2.6% 

5 

10 

+4.0%, -2.2% 

5 

11 

+3.8%, -2.6% 

5 

Average 
w/o Run 5 

+3.6%, -2.8% 
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Table 3 


Location of Test Data for Each Test Run 



Run 

Run 

Run 

Run 

Run 

Run 

Run 

Run 

Run 

Run 

Run 


1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

ii 

Test Conditions 

T4 





















and Mass Flows 

T5 

T5 

T5 

- 

T5 

T5 

T5 

- 

T5 

T5 

- 


Fll 

Fll 

Fll 

- 

Fll 

F12 

F12 

- 

F12 

F12 

- 

Static Pressures 

F13 

F13 

F13 

T6 

F13 

F14 

F14 

T6 

F14 

F14 

T6 

Centerline Veloci- 

F15 

F15 

F15 

_ 

F15 

F16 

F16 

- 

F16 

F16 

- 

ties and Temper- 












atures 












Velocity Profiles 

F17 

F18 

F19a 


F20 

F21 

F22 

_ 

F23a 

F24 

- 

1 



F19b 






F23b 



Temperature 












Profiles 

- 

- 

F25 


- 

- 

- 

- 

F26 


— 




F27 



F27 






Eddy Viscosity 

- 

- 

F28 

- 

- 

F28 

- 

- 

- 

- 

“ 

Sensitivity 



F29 



F29 






Flow Rate 



F37 



F27 




. 

_ 

Sensitivity 













T Stands for Table 
F Stands for Figure 
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Summary of Experimental Test Conditions and Flow Rates 


Mass- 

Momentum 

Stagnation 

Pressure 
at x = 10. 5” 
Dsia 

16.15 

16.19 
15.91 

15.86 

16.23 

16.20 

16.06 

15.96 

16.05 

Mixing 

Section 

Throat 

Size 

1. 25" 
1. 875" 

Flow 

0 

1 

t- 

. ^ COtOCOCQooC^ODWt-OC-’ 

» 

£ . , . 

Secondary 

Flow 

a 

S o 
2 8 
5 

m OWNOOiclOOiiONt-aO 

> cocaoicocoeocoeo ^^^ 1 

Mixing 
Section 
Flow 
lb/sec in, 

l 

P CO rH tr- <0 Tf 00 Tf iH ^5 to 

> *P w 

odd d d d d d d Ct a 

0) 

N 

N 

0 0) 

**a 

01 P$ 

.2 

i 

S ! 

W N 

0, 0780 
0.0782 
0. 0750 
0. 0874 
Q. 0787 
0. 0882 
0. 0884 
0. 0782 
0. 0884 
0. 0874 
0. 0876 

00 O 
2 U-* 

s 

lb/sec 

0. 6240 
0. 6263 
0. 5997 
0. 6993 
0. 6293 
0,7056 
0.7069 
0.6256 
0. 7072 
0. 6993 
0.7006 

Atmo- 

spheric 

a 

Stf 

©o 

— i 

T 

a 

538 

543 
553 
553 

544 
547 
543 
543 
550 

547 

548 

. 

Baro- 

metric 

Pressure 

psia 

ChOrHr-CrHOOOCOOOD 
tDtDtOtOC— OOOOOOt— t*-t£> 

fi tt*B9l««SSS 

Cl 

t— It— 1?— lr-4Wi— It— It-Ht-^t— (t— ( 

Nozzle 

Throat 

Co- 

efficient 

Tf m 

- t-^ooaoo^^coLOOco 

M tOt£>t-t-aOt>t>tDI>t-t- 

O OJOiOlOOiOOOOOid 

• ..**«••••• 
ooooooooooo 

£, t5 

S3 § 

N U >-i ^ 

0.9688 

g § 
|H° 

T 

N 

641 
637 
706 

659 

648 

649 
647 

642 
644 

660 
652 

Nozzle 

Pressure 

psia 

_ CnO»HWr-fOOOCOOCp 

^ tococotot-aoooooe-c-to 

Cl * .««••••*!• 

eocoeQCoeoeococoaocQCO 


Si 

PS Z 

Symbol 

1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 
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Table 5 


Comparison of Experimental and 
Analytical Mass Flow Rates 


Run 

No. 

Throat 

Width 

(inches) 

Mixing 
Section 
Mass Flow 
Rate From 
Traverse 
Data 

lb/ sec. in. 

CD 

Mixing 
Section 
Mass Flow 
Rate From 
Orifice 
Data 

lb/ sec. in. 
© 

Percent 

Difference 

In 

Measured 

Data 

© - © 

Analytical 
Mass Flow 
For 

Best Static 
Pressure 
Match 
lb/ sec. in. 
© 

Comparison 

Of 

Traverse 

To 

Analytical 
Mass Flow 

CD - © 

Comparison 

Of 

Orifice 

To 

Analytical 
Mass Flow 

© - © 

HSH 

© 

HgHI 

1 

1 


0.408 


.. 

0.382 

+ 6.4% 


2 

■ 


0.341 

0.320 

+ 6.7% 

0.322 

+ 5.6% 

- .6% 

3 

l! 


0.357 

0.3535 

+ 1.0% 

0.351 

+ 1. 7% 

+ .7% 

5 

1 


0.384* 

— 

— 

0.395 

- 2.9% 

— 

6 

I 


0.434 

0.417 

+ 4.0% 

0.420 

+ 3.2% 

- .7% 

D 

1 


0.458 

| 

0.447 

| 

+ 2.5% 

0.443 

+ 3.3% 

+ .9% 

H 



0.501 

I 

1 

— 

0.485 

+ 3.2% 

— 

10 

1 


0. 525 

Hi 

— 

0.508 

+ 3.2% 

HI 


* Transducer Battery May Have Been Going Bad During This Test 
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Table 6 



Tabulation of Static Pressures For 
Runs 4, 8, and 11 


Distance 
From Nozzle 

Run No. 

4 

8 

11 

Discharge 

Nozzle Pressure psia 

35 , 61 

31.80 

35.68 

inches 

Throat Height 

1 . 25 " 

1 . 875 " 

1 . 875 " 


Wall Static Pressure in Inches of Water Gage 

- 1 . 62 


- 8.05 

- 6.8 

- 10.4 

- 1.25 


- 11.2 

- 9 . 1 

- 13.7 

- 0.87 


- 13.7 

- 10.7 

- 16.3 

- 0.46 


- 17.1 

- 13.0 

- 19.5 

+ 0.03 


- 19.5 

- 14.2 

- 21.5 

+ 0.56 


- 19.5 

- 13.9 

- 20.9 

+ 0.99 


- 20.7 

- 14.2 

- 21 . 5 

+ 1.50 


- 21.5 

- 13 . 9 

- 22.4 

+ 2.00 


- 28.0 

- 16.5 

- 26.6 

+ 2.50 


- 32.2 

- 17.7 

- 28 . 6 

+ 3.00 


- 35.7 

- 18.6 

- 30 . 1 

+ 3.50 


- 37.5 

- 18.3 

i 

to 

o 

h -* 

+ 4.00 


- 39.5 

- 18.6 

- 31.0 

+ 4.50 


- 40.7 

- 18.3 

- 31.3 

+ 5.00 


- 41.9 

- 18.3 

- 31 . 6 

+ 5 . 50 


- 42.7 

- 18.0 

- 31.9 

+ 6.00 


- 43. 1 

- 18.0 

- 32 . 5 

+ 7.00 


- 48.7 

- 18.6 

- 33.6 

+ 8 , 00 


- 54.3 

- 19.2 

- 35.4 

+ 9 . 00 

; 

- 55.8 

- 18.5 

- 35.4 

+ 10.50 


- 58 . 1 

- 17.7 

- 35 . 1 

+ 12.00 


- 42.2 

- 13.0 

- 28.6 

+ 13.00 


- 30. 1 

- 9. 1 

- 23.3 

+ 15.00 


- 13 . 2 

- 2.5 

- 15 . 1 

+ 17.00 


- 2.3 

+ 2.2 

- 8.3 

+ 19.00 


+ 4.9 

+ 5.9 

- 5 . 1 

+ 21.00 


+ 9.9 

+ 8 . 6 

+ 0.7 

+ 23.00 


+ 13.2 

+ 10.9 

+ 3.7 
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Figure 1 

Assembly Sketch of Two Dimensional Ejector Test Rig 







Figure 2 

Picture of Primary Nozzle 
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Figure 3 

Picture of Nozzle Positioned in the Mixing Section 
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1 



E)© © 0 © © (f 





Bleed 

Valve 
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Schematic of Experimental Layout 
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Figure 7 

Picture of Right Side of Ejector Rig 



Figure 8 

Picture of Left Side of Ejector Rig 
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Mixing Section Static Pressures 










5 


0.5 


Figure 12 


W m - Mixing Section Total Mass Flow - lb/sec. in. 

Comparison of Experimental and Analytical Mass Flow Rates for 
Runs 6, 7, 9, and 10 




Wall Static Pressure - Inches of Water Gage 










Wall Static Pressure - Inches of Water Gage 


Solid Lines are Analytical Results 


+20 



Figure 14 Wall Static Pressure Distributions for Mixing Section 
with 1. 875" Throat 
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Maximum Velocities - fps 



200 ^ 
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50 
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Figure 16 Maximum Velocities for 1.875" Throat Mixing Section 

P N = 21. 0 psig 
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Centerline Temperature 






Figure 18 Velocity Profiles for Run 2 for 1.25" Throat Mixing Section 
P N * 17. 0 psig T n = 177° F, W N * 0.0782 lb/ sec. in. 
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Velocity - fps 



Top 

Wall 


Fractional Distance From Centerline 


Bottom 

Wall 


Figure 19b Velocity Profiles for Run 3 for 1. 25” Throat Mixing Section 
p N « 17.0 psig, T n * 246* F, W N * 0.075 lb/ sec. in. 
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Figure 20 Velocity Profiles for Run 5 for 1.25” Throat Mixing Section 
P N - 17. 0 psig, T n * 188° F, W n « .0787 lb/ sec. in. 







Velocity - fpe 



Figure 21 Velocity Profiles for Run 6 for 1.875” Throat Mixing Section 
^N*21. 0 psig, Tj^^'F, W N » .0882 lb/sec. in. 
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Velocity - fps 


1300 


Analytical W = 0.443 lb/sec. in. 



Figure 22 Velocity Profiles for Run 7 for 1.875” Throat Mixing Section 
P N = 21. 0 psig, T n = 187°F, W N = . 0884 lb/sec. in. 
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Analytical W M = 0.485 lb/sec, in. 



Figure 23a Velocity Profiles for Run 9 for 1.875" Width Mixing Section 
= 21. 0 psig, = 184°F, W.^. = . 0884 lb/sec. in. 
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Wall 


Figure 23b Velocity Profiles for Run 9 for 1.875 fl Throat Mixing Section 
P M = 21. o psig, T xt = 184°F, W XT - . 0884 lb/sec. in. 







Velocity - 
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Wall 


Fractional Distance From Centerline 
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Wall 


Figure 24 Velocity Profiles for Run 10 for 1.875" Throat Mixing Section 
P N = 21. Opsig, T^ = 200°F, W^ = .0874 lb/sec. in. 
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©□ - Test Data 



Wall Fractional Distance from Centerline Wall 

Figure 26 Temperature Profiles for Run 9 for 1.875" Throat Mixing Section 


125 




Constants Eq. 20 is 0.094, Eq. 21 is 0.120 
Constants Eq. 20 is 0.08, Eq. 21 is 0.108 
Constants Eq. 20 is 0.094, Eq. 21 is 0.120 


lb/ sec. in. 


Hun 6 


lb/ sec. in. 


\> 




lb/ sec. in 


Axial Distance From Nozzle Discharge - Inches 

Figure 27 Wall Static Pressure Sensitivity to Flow Rate and 
Eddy Viscosity for Run 3 and Run 6 








Constants - Eq. 20 is 0. 08, Eq. 21 is 0. 108 

Constants - Eq. 20 is 0. 094, Eq. 21 is 0. 120 



Figure 28 Centerline Velocity and Temperature Sensitivity to Eddy Viscosity 
for Run 3 and Run 6 
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Centerline Temperature 










Throat Wall Static Pressure - Inches of Water Gage 






